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ABSTRACT 


This  study  develops  correction  factors  for  currently  used  stress 
intensity  factor  equations  to  more  accurately  predict  stress  intensity 
factors  for  a  corner-crack  emanating  from  a  hole  as  it  transitions  to 
a  uniform  through-the-thickness  crack.  These  correction  factors 
resulted  in  an  increase  in  accuracy  for  total  life  prediction  and  much 
better  correlation  between  analytical  stress  intensity  factor 
predictions  and  ex  per  imenta  1  results  in  the  transit  ion  region  for 
_!< )L24— T35 1  and  7075-1651  aluminum  alloys.  An  experimental  program  was 
undertaken  to  generate  all  2024-1351  aluminum  Lest  data  used  in  this 
investigation.  The  7075-T651  aluminum  test  data  was  generated  during 
earlier  work.  Correction  factors  developed  by  Opel  for  7075-T651 
aluminum  were  evaluated  for  2024-T35I  aluminum  and  found  to  be 
excess! vely  conservat ive. 

ilartranft  and  Sih  suggested  the  stress  intensity  factors  be 
evaluated  at  an  imaginary  surface  at  an  angle  0  a wav  from  the  front 
and  hole-bore  surfaces.  These  stress  intensity  factors  could  then  be 
used  in  the  life  prediction  models  for  Life  from  a  corner-crack  until 
hack  surface  penetration.  Evaluating  stress  intensity  factors  at  0= 
Inland  S0°  eliminate  surface  boundary  effects  caused  by  fabrication 
processes  like  cold  rolling  and  hole  drilling,  thereby  improving  life 
predict  ions  to  back  surface'  penetration  significantly.  It  is  shown 
that  with  ‘a  =  10 J  and  TO  ’  ,  no  other  corrections  are  necessary  for  life 
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prediction  tn  back  surface  penetration. 

A  transition  region  is  postulated  to  effect  life  prediction  from 
back  surface  penetration  until  final  fracture.  Correlations  were  made 
be  t  ween  experimental  results  tor  2024- mi  and  7075-T651  aluminum  test 
data  and  the  Crunch  1  i  near  i /.at  i  on  ot  the  Bowie  solution  for  a 
t norouyh-! he-t h ic kness  crack.  These  correlations  were  plotted  from 
!),lt  ''  sur:  ”  p  iwuu't  rat  ion  until  final  fracture.  The  plot  yields  the 
,>n,!  r,’H'  ,r  iftsition  region  and  the  required  correction  factors.  The 

;i:,‘  aroint  ion  model  incorporating  these  transition  correction 

'  MO  siiown  to  be  the  most  accurate  and  versatile  ot  all  models 

mo.  -  St  i  eat  ed  . 

i  I predictions  were  made  using  the  new  model,  an  instantaneous 
-odol  (which  assumes  a  through-crack  at  hack  surface  penetration), 
OpcC’s  model,  for us sat ' s  model  and  the  Collipr i est-Ehret  model. 
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1.0  INTRODUCTION 


1.1.  Mot  i  vat  ion 

Cracks  Ln  aerospace  structural  components  continue  to  cause 
problems  ranging  from  premature  retirement  of  parts  to  catastrophic 
fa  i  Lures.  These  problems  result  in  the  loss  of  millions  of  dollars 
annual lv  to  the  aerospace  industry  and  have  caused  loss  of  life  in 
certain  cases.  A  good  deal  of  progress  has  been  made  over  the  past 
several  years  to  eliminate  the  conditions  leading  to  crack  initiation 
and  growth.  However,  flaws  due  to  materiaL  defects,  manufacturing 
methods,  and  in-service  conditions  will  always  be  present.  Therefore, 
in  order  to  correctly  predict  component  fatique  life,  a  life- 
prediction  method  which  assumes  the  presence  of  flaws  at  the  onset  of 
the  component  operational  usage  is  required. 

In  1071  ,  the  1'nited  States  Air  Force  (USAF)  sponsored  an  in-depth 
review  ol  aerospace  structural  failures  (Gran  et  al, [1]).  This  review 
showed  that  the  predominate  failure  mechanism  was  cracks  emanating 
from  fastener  holes.  In  fact,  over  one-third  of  all  failures  studied 
were  due  to  cracks  emanating  from  fastener  hoi  .  Ac  a  result  ol  this 
review  the  I'SAF  recently  adopted  a  damage-tolerant  design  philosophy 
t<>  ensure  that  catastrophic  failures  are  not  caused  bv  cracks 
emanat  i ng  from  fastener  holes.  This  philosophy  was  introduced  by  Wood 
and  Fugle  1 2 1  and  is  based  on  Linear  'clastic  Fracture 


Median  ics  ( LKl’M ) . 


This  fracture  mechanics  analysis  of  cracks  at  fastener  holes 
requires  knowledge  of  the  stress  intensity  factor,  K,  for  all  crack 
geometries  of  interest.  Since  the  crack  length  within  the  bore  of  the 
hole  cannot  be  accurately  measured,  a  combination  of  analytical  crack 
solutions  must  be  used.  Exact  analytical  expressions  for  stress 
intensity  factors  have  been  derived  for  cases  involving  very  small 
flaws  in  infinite  (very  large)  elastic  bodies  (e.g.,  Williams  [3], 
[4]).  These  solutions  have  to  be  modified  for  finite  problems  by  the 
inclusion  of  correction  factors.  These  correction  factors  are 
generally  derived  using  experimental  data  and  backtracking  techniques, 
finite  element  approximations,  or  any  of  various  other  approximate 
methods  (see  e.g.  see  references  [5]-  [11],  which  are  used  in  this 
work).  The  stress  intensity  factors  are  functions  of  the  test- 
specimen  geometry,  loading,  and  crack  geometry.  Many  works  published 
over  the  past  several  years  have  dealt  with  estimating  and/or 
i ulculat  ing  stress  intensity  factors,  and  the  most  significant  of 
these  works  are  discussed  in  the  next  section. 

1 . 2  Background 

1.2.1.  Previous  Work  on  Stress  Intensity  Factors  for  Corner- 
Cracks  Emanating  from  a  Hole.  In  1956,  Bowie  [12]  made  the  first 
major  contribution  to  the  problem  of  cracks  emanating  from  holes, 
specifically,  cracks  emanating  from  open  holes  in  plates.  Using 
complex  variable  methods,  Bowie  solved  the  two  dimensional  problem  of 


single  .uni  double  through-cracks  at  an  open  hole.  Bowie's  work  is 
signil leant  in  that  although  it  is  not  an  exact  solution,  it  is 
universally  used  as  a  basis  for  comparison  with  other  studies  of 
cracks  emanating  from  holes.  More  importantly,  Bowie's  solution  is 
also  widely  used  to  establish  the  accuracy  of  other  methods  which 
attempt  to  solve  more  complicated  crack  problems. 

Bowie's  solution  was  improved  upon  in  1972  bv  Liu  [13].  He 
studied  the  quar t er-el 1 i pt ica 1  crack  at  a  hole  in  a  plate  and 

.  i  imposed  a  number  of  solutions  for  his  approximation.  He  used  the 
Smith  et  a  1 .  solution  [14]  to  account  for  the  hole  surface  and  the 
t rant  surface  of  the  specimen  and  Kobayashi's  solution  [15]  as  a 
correction  tor  the  back  surface.  These  approximations  were 
superimposed  on  Bowie's  [12]  two  dimensional  approximation  of  the 
three-dimensional  hole  effect.  This  did  not,  however,  result  in  a 
solution  and  only  gave  stress  intensity  estimates  at  a  point  on  the 
periphery  of  the  crack,  midway  between  the  front  of  the  specimen  and 
the  surface  of  the  hole. 

A  unique  approach  to  the  solution  of  the  single  through-crack 
problem  was  presented  by  Tweed  and  Rooke  in  1973  [16].  They  derived 
st ress- i n tens i t v  factor  relations  using  a  Mellin  integral  transform 
technique.  This  method  provided  a  substantial  improvement  over  the 
Bowie  solution,  especially  for  small  crack  lengths. 


Grandt  |  10]  used  the  principle  of  superposition  to  develop  a 


stress  intensity  factor  solution  for  large  plates  containing  radial 
holes  with  through-the-t h ickness  cracks  loaded  with  arbitrary  crack- 
face  pressure.  Grandt's  method  is  based  on  work  by  Rice  [17],  who 
showed  that  once  the  displacement  field  and  stress  intensity  factor 
are  known  for  one  geometry  and  loading,  Kj  the  stress  intensity  factor 
for  mode  I  deformation  may  be  obtained  for  any  other  symmetric  loading 
applied  to  the  same  crack  geometry.  Grandt  [10]  used  a  remote, 
uniform  tensile  stress  as  the  known  loading,  crack-tip  stresses 
obtained  by  Bowie  [12],  and  stress  intensity  factor  results  bv  Paris 
and  Sih  [Id]  to  develop  his  solution  tor  through-the-t h ickness  tracks 
at  holes.  Grandt  developed  solutions  for  t  number  ol  different 
fastener  conf i gurat ions. 

In  1078,  Kullgren  et  al.  [10]  extended  and  refined  the  numerical 
solution  introduced  by  Browning  and  Smith  [20]  that  was  based  on  a 
modification  of  the  Sell  wart  /-Neumann  alternating  met  hod  lor  oh'  lining 
st  ress-i  ntensi  t  v  factors  for  a  sei:ii-<  ircular  <  r.u  k  it  i  hole. 

k’uligren's  technique,  labeled  the  finite  element  1 1 1  et  nat  i  ug  . 1'iog, 

is  based  on  finite  element  .tpproxin.it  ions.  Kullgren  .  •  .nsi  del  eg 
quarter-elliptical  cracks  emanating  t  run  open  holes  i  :  i  plate-., 
ktillgren  and  Smith  |2i|  extended  this  method  to  single  sen  i  -e !  !  . -a  < .  t; 
mil  double  quar  t  er-e  1  1  i  pt  i  c.a  1  c  r.u  ks  en.iii.it  i  n  >  '  p  «n  t  i  si  ener  gele-.. 
fills  work  iias  demons!  rat  ed  th.it  the  linite  element  i  1 1  ern.tt  :  ng  ".et'.'g 
is  versatile.  It  mii  be  adopted  tor  m.inv  liMereiu  tomplex 


t*i  jint*  t  r  i  us ,  load  ini;  conditions  and  crack  orientations.  The  only 


drawback  of  this  method  seems  to  be  the  difficulty  involved  in 
modi! vim;  the  computer  codes. 

A  rather  significant  contribution  to  the  area  of  three- 
dimensional  cracks  subjected  to  uniform  tension  was  made  by  Newman  and 
Ra j u  I1)]  in  1  ‘ >H 1 .  They  studied  several  crack  configurations  and 
considered  several  parameters  such  as  crack  depth,  crack  length, 
specimen  thickness  and  hole  radius.  Thev  considered  solutions  due  to 
Howland  [  J2  J  for  stress  concentration  at  a  hole  and  Iada  et  al.  |  _’ 5  ] 
for  crack  eccentricity  corrections  as  a  basis  for  their  finite  width 
correction  factors.  Newman  and  Raju  also  evaluated  Shah's  (onvorsion 
factors  that  were  applied  to  approximate  double-crack  solutions  us  me 
si  nt>l  e-crack  results.  These  factors,  which  were  derived  In  usi-- 
Green's  function  to  find  stress  intensity  factors  at  a  hole  : •  ,  .  • 

under  uniaxial  tension,  were  found  to  be  in  nood  agreement  wit  h  '  > 
recall  ts  of  Kulluren  and  Smith  [21  |. 

Solutions  tor  analytical  stress  intensity  factors  t<a  ■  .a 
cracks  at  ■>  hole  due  to  Newman  and  Raju  | 6  J  ,  Shall  |  24  |  and  . 
were  evaluated  by  Heckel  and  Rudd  [2a]  usi  nu  correlation1, 
experimental  stress  intensity  factors  obtained  nsiii"  the  I  !.:■<■  - 
Anderson  bat  kt  rack  i  n<>  technique  |  26  ] .  Thev  found  that  Newman  mu 
Ra ju  s  solution  provided  excellent  analytical  'experiment.:! 
correlations  ot  the  stress  intensity  factor,  crack  crowtl,  i  it  ev  mu 


'  ra<  k  shape  <  liarnies.  Shah's  solut  ion  produced  e\c  e!  I  f  ■  n  t 


•  t  h'ss 


intensitv  factor  correlations  and  life  predictions.  I'hey  also 


detemined  that  both  of  Liu's  one-dimensional  and  two-dimensional 
solutions  vielded  unconser vat i ve  stress  intensity  factors  and  crack 
growth  rates. 

In  1983,  Schi jve  [27]  used  data  collected  by  Raj u  and 
Newman  (28,  29]  for  semi-  and  quarter-el  1 i pt ical  corner  cracks  at 
holes  to  develop  interpolation  techniques  for  cracks  with  other 
dimensions,  lie  noted  that  real  cracks  will  not  always  have  elliptical 
shapes,  hut  as  a  first  approximation  t  lie  semi-  and  quarter-elliptical 
r  t*  1  trout  issumpt  ion  is  a  valid  one.  However,  Schi  jve  noted  that  in 
•  1 1  ,  'he  ic  met  r  i  <  ratios,  a/c,  a/t,  R/t  (where  a  is  crack  depth, 

i  ■  1  eiiet  h ,  »  is  spe<  t  men  thickness  and  R  is  hole  radius)  will 

•  r  :•  » tiose  used  bv  Ra  ju  and  Newman .  Schi  jve  remedied  this  by 
it  c  interpolation  methods  judiciously  chosen  for  qiven  crack 
•  .  Lie  results  were  compared  to  the  Raju-N’ewman  solutions 

c  '  ■  ciec  to  within  22 . 

.  .  r  in-. i  i  ion  .a  i  Corner  Track  at  a  Hole  to  a  Through- the- 

t  i-  •  ■  he  transition  res’ ion  starts  when  the  corner  crack 

.  ■  o  t  the  hole  nearly  reaches  the  point  of  hack 

>  i*  ■  ui  1 1 1 . i  ends  when  the  crack  becomes  a  uniform  through- 
-■  r  i.'.  'Determining  the  boundaries  of  this  rest  ion  is  a 

i  '  !■'  in'  .t  i  nis  work  and  will  he  discussed  in  detail  in  a 

i  i  i  r  e n  t  I  v ,  procedures  for  dealing  with  this 
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correction  lac  tors  to  compensate  lor  transition  effects  to  totally 
ignoring  the  transition  region.  The  American  Society  for  Testing  and 
Materials  (ASTM)  is  trying  to  remedy  this  rather  inconsistent  approach 
by  sponsoring  "round-robin"  competitions  for  predicting  the  life  of 
specimens  to  determine  the  accuracy  of  current  procedures.  these 
round-robin  competitions  have  aided  greatly  in  the  rapid  advances  in 
the  development  of  transition  region  criteria.  ASTM  published  some  of 
these  results  dealing  strictly  with  part-through  crack  fatigue  life 
predictions  in  1979,  (see  Chang  [30]). 

Peterson  and  Vroman  [31]  developed  a  computational  scheme  to 
determine  stress  intensity  factors  for  initial  part-through  crack 
defects.  Instead  of  using  the  number  of  loading  cycles  as  the 
independent  variable,  as  in  most  current  solutions,  their  computer 
code  uses  very  smaLl  increments  of  crack  si/.e  to  calculate  the  number 
oi  loading  cycles  required  for  each  increment  of  crack  growth.  The 
code  uses  equations  developed  by  Forman  [32]  to  calculate  stress 
intensity  factors  at  the  surface  and  along  the  crack  front.  K,  the 
stress  intensity  factor,  is  determined  as  a  I  unction  of  crack  depth  to 
specimen  thickness  ratio,  crack  aspect  ratio,  crack  depth,  and  crack 
half-length.  Crack  growth  rates  were  calculated  bv  substituting  the 
stress  intensity  factors  into  a  Paris-type  equation  (33],  When 
Peterson  and  Vroman  first  presented  their  approach  in  l‘>72,  it  was 
found  to  be  more  conservative  than  the  Kohavash i -Moss  solnt  ion  |  13]. 
Further  invest i gat  ion  in  1974  concluded  that  the  Kohavash i -Moss 


solution  was  more  accurate.  However,  it  was  also  concluded  that  the 
Peterson-Vroman  procedure  was  attractive  due  to  its  ease  of 
application.  A  1476  modification  to  the  Peterson-Vroman  approach 
provided  a  good  approximation,  and  made  this  method  even  more 
at  tract ive. 

The  common  procedure  for  dealing  with  the  transition  of  a  corner 
crack  at  a  hole  to  a  through-crack  is  to  assume  a  corner  crack  until 
back-surface  penetration  and  then  assume  a  through  crack  until 
fracture.  This  procedure  is  quite  conservative.  Johnson  [34] 
compensated  for  this  by  suggesting  that  the  crack  be  considered  a 
t hrough-the-thickness  crack  only  after  the  back-surface  crack  length 
is  ()0-,  of  the  front-surface  crack  length.  Johnson  based  his  idea  on 
an  imaginary  crack  depth  found  by  letting  the  crack  grow  in  the  same 
elliptical  shape  it  had  when  it  penetrated  the  back  surface.  This 
assumption  seems  to  be  in  good  agreement  with  the  trends  of 
experimental  data  [5,  6). 

Brussat  and  Chiu  [8]  considered  the  crack  in  the  transition 
region  lo  be  a  combination  of  a  surface  flaw  and  a  through-the- 
Lhickness  crack.  Their  criterion  was  based  on  crack  depth,  length, 
and  specimen  thickness  for  a  quarter-circular  corner  flaw  at  a  hole. 
Brussat:  and  Chiu  developed  transition  region  correction  factors  based 


mi  engineering  judgment  rather  than  on  actual  mathematical  results. 


Tlie  most  commonly  used  procedure  for  predicting  the  Life  of  a 
specimen  considering  the  transition  region  of  a  corner-crack-at-a-hole 
is  the  one  discussed  by  Engle  in  [11].  This  procedure  has  been 
adopted  by  the  USAF  as  the  current  damage-toierant  philosophy  and 
forms  the  basis  and  starting  point  for  this  thesis.  It  uses  Newman 
anti  Raju's  corner  crack  solution  until  the  crack  depth  is  equal  to  the 
plate  thickness,  and  then  uses  Grandt's  [10]  linearization  of  Bowie's 
solution  [12]  for  a  through-the-thickness  crack  near  a  hole.  However, 
this  apnroach  ignores  the  portion  of  the  transition  region  extending 
from  back-surface  penetration  to  where  the  crack  becomes  a  through- 
the-thickness  one.  The  results  are  therefore  conservative. 

In  ll)<33,  Opel  [5]  suggested  modifications  to  the  method  discussed 
by  Engle.  Opel  improved  the  total  life  predictions  by  approximately 
IS’7,.  His  work  also  resulted  in  a  better  correlation  between  the 
analytical  stress  intensity  factor  prediction  and  the  experimental 
results  by  developing  transition-region  correction  factors  for  corner- 
t.  racks-at-a-hole  tests.  Opel's  work  was  performed  in  three  phases. 
First,  correlations  were  made  between  experimental  results  obtained  by 
Grandt  and  Snow  [35]  for  Polymethylmethacrylate  (P1IMA)  testing  and  the 
Newman-Ra ju  [0]  three-dimensional  stress  intensity  factor  equation  for 
a  single  corner  crack  at  a  hole.  These  correlations  became 
unconservative  when  the  normalized  crack  depth  (a/t)  reached  a  value 
ot  0.75,  and  this  was  defined  as  tiie  beginning  of  the  transition 
region.  Next,  correlations  were  made  between  experimental  results 


obtained  by  Heckel  and  Rudd  125,  36]  for  7075-T651  aluminum  testing 
and  the  Grandt  linearisation  of  the  Bowie  solution  for  a  through  crack 


emanating  from  a  hole.  These  correlations  were  found  to  be 
unconservative  until  the  normalised  crack  depth  (c/R)  reached  a  value 
of  2.5;  beyond  that  point  the  correlations  became  conservative.  This 
point  was  defined  as  the  end  of  the  transition  region.  Finally, 
stress  intensity  factors  for  the  transition  region  were  developed. 
Along  the  bore  of  the  hole,  a  second-order  polynomial  regression 
curve-fit  was  used  for  the  normalized  crack  depths  (a/t)  found  in  the 
transition  region.  The  resulting  correction  factor  was  multiplied  by 
the  stress  intensity  factors  found  using  the  Mewman-Ra ju  solution.  A 
first  order  polynomial  regression  curve-fit  was  found  to  give  accurate 
correction  factors  for  the  N'ewman-Ra ju  stress  intensity  factors  along 
the  front  surface.  A  second-order  polynomial  regression  curve-fit  was 
used  to  find  the  transition  correction  factors  for  the  Grandt 
1 i near izat i on  of  the  Bowie  solution  along  the  back  surface  until  final 
fracture.  Opel  proceeded  to  make  life  predictions  using  his  new 
model,  the  model  referenced  by  Engle  [11],  and  the  Brussat  model  [3|. 
Opet's  model  proved  to  be  the  most  accurate. 

1 . 3  Scope 

An  experimental  program  was  undertaken  to  investigate  the  crack- 
growth  behavior  in  rorner-r rack-at-a-hole  specimens  of  2024-1351 
aluminum  alloy.  01  particular  interest  was  an  assessment  of  the 
effectiveness  of  a  procedure  introduced  bv  Opel  [5]  for  predicting 


rrack  growth  rates  in  similar  specimens  of  7075-T651  aluminum  alloy  to 
predict  crack  growth  in  other  materials.  He  used  experimental  data 
obtained  from  PUMA  and  7075-T651  specimens  to  derive  correction 
iactors  !'•  -  existing  solutions  which  were  shown  to  yield  better 
predictions  for  these  materials.  Just  as  in  Opel's  work,  a  part- 
elliptical  corner  crack  at  a  hole  is  considered.  The  corner  crack 
center  is  located  at  the  intersection  of  the  hole  wall  and  the  front 
surface.  The  crack  lies  in  a  plane  perpendicular  to  the  axis  of 
loading.  The  finite  plate  geometry  and  loading  conditions  are  such 
that  the  plane  of  the  crack  is  a  plane  of  symmetry  for  the  open  hole 
problem.  The  loading  is  a  remote  uniaxial  constant  amplitude  loading, 
and  the  initial  crack  eccentricity,  a/c ,  is  greater  than  one. 

The  specific  objectives  of  this  thesis  are  to: 

1.  Evaluate  the  applicability  of  Opel's  [5]  correction  factors 
to  2024-T351  aluminum,  and  should  they  prove  to  yield 
inaccurate  predictions,  calculate  new  correction  factors  for 
the  stress-intensity-factor  models. 

2.  Evaluate  the  increase  in  accuracy  suggested  by  ilartranft  and 
Sih  137]  for  life  predictions  if  the  stress  intensity  factors 
are  calculated  at  an  angle  0  from  the  surfaces  as  indicated 
in  Fig.  I  . 

3.  Determine  if  transition  correction  factors  for  the  region 
from  corner-crack  until  back  surface  penetration  are  needed 
it  the  ilartranft  and  Sih  suggestion  is  used  as  mentioned  in 


4.  Detenu  no  the  transition  region  and  new  transition  correction 
factors  for  the  development  of  an  improved  life  prediction 
mode  1 . 

5.  Evaluate  the  ability  of  existing  stress  intensity  factor 
models  to  predict  the  behavior  of  a  corner-crack-at -a-iiole  in 
2024-T351  aluminum  as  it  transitions  into  a  uniform  through- 

t he-thickness  crack  (Fig.  1). 

(>.  Assess  the  accurarv  of  the  corrected  model  by  comparing  the 
total  life  predictions  obtained  using  these  models  to  results 
of  other  life-predictions. 


Stress  intensity  factor  correlations  were  performed  for 
experimental  results  obtained  at  the  Air  Force  Flight  Dynamics 
Laboratory  for  2024-1351  and  7075-T651  aluminum,  constant-amplitude¬ 
loading  testing  and  the  Orandt-Bowie  [10,  11,  12]  stress  intensity 
solution  tj r  a  t hrough-t he-t h i ckness  crack.  Opel's  model  has  not 
proved  to  be  adequate  for  a  new,  more  ductile  material  (2024-T151 
aluminum).  Examination  of  Table  1  reveals  that  although  the  Opel 
model  shows  a  9%  increase  in  accuracy  over  instantaneous  model  for 
2024-T35I  aluminum,  the  total  life  predictions  arc  still  300 
conservative.  Therefore,  a  more  accurate  model  will  be  investigated. 

Those  correlations  were  [(lotted  from  back-surface  penetration  to 
final  fracture.  Correlations  were  performed  for  both  front  and  back 
surfaces.  These  plots  were  used  to  develop  transition  correction 


inn  i  mini 


factors  for  the  region  from  back  surface  penetration  to  the  end  of  the 
transition  region,  i.e.,  to  the  point  where  a  uniform  through-the- 
tliickness  crack  deveLops.  The  end  of  the  transition  region  was 
identified  at  the  point  where  the  back  surface  crack  length  is  nine- 
tenths  of  the  front  surface  crack  length  (i.e.,  c'  =  .9c,  as  in 
Johnson [ 34 ) )  (More  details  on  the  definition  of  the  transition  region 
will  be  given  in  later  sections.)  These  correlations  and  transition 
correction  factors  are  compared  to  those  developed  and  presented  Dv 
Opel  [  5 J  and  Opel,  Rudd,  and  Ilaritos  [6). 

Finally,  total  life  predictions  are  made  from  a  specified 
initial  crack  length  to  final  fracture  using:  a)  a  transition 
criterion  developed  in  this  paper,  b)  an  instantaneous  transition 
criterion  presented  in  Chang  [20] 

c)  at  the  transition  criteria  developed  by  Opel  [5],  d)  the  criterion 
developed  by  Brussat  et  al  [8],  and  e)  the  corrections  developed  by 
Collipriest  and  Ehret  [7].  The  criterion  developed  in  this  paper  is 
then  evaluated  by  comparing  its  predictions  to  those  of  Chang,  Brussat 
et  al.,  Opel  and  Collipriest  and  Ehret,  as  well  as  to  the  experimental 
results  obtained  in  the  course  of  this  investigation. 


2.0  EXPERIMENTAL  PROGRAM 


The  entire  experimental  program  was  conducted  at  the  Air  Force 
Wright  Aeronautical  Laboraties,  Flight  Dynamics  Laboratory,  Wright- 
Patterson  Air  Force  Base,  Ohio.  Two  distinct  sets  of  data  were 
generated  for  this  study:  mater ial ' character ization  data  and  corner- 
crack-at-a-hole  data,  both  for  2024-T351  aluminum.  The  material- 
characterization  tests  included  both  crack-growth-rate  and  fracture- 
toughness  tests.  The  James-Anderson  backtracking  technique  was 
applied  to  these  data  in  order  to  obtain  experimental  stress  intensity 
factors.  These  factors  were  used  in  the  life  predictions  for  the 
corner-crack-at-a-hole  tests.  The  corner-crack-at-a-hole  test  results 
were  used  to  develop  a  new  transition  criterion  and  evaluate  the 
accuracy  of  this  new  criterion  by  comparing  it  to  the  results 
reported.  The  details  of  all  tests  conducted  and  the  procedures 
followed  are  described  in  the  next  sections. 

2 . 1  Material-Characterization  Tests 

Material-characterization  tests  were  conducted  using  center- 
cracked  tension  specimens.  (Specimen  dimensions  are  given  in  Fig.  2) 
The  crack-growth-rate  tests  used  three  sets  of  duplicate  specimens, 
with  maximum  stress  a  max  -  10.0  ksi  and  stress  ratio  R  =  0.1,  0.5, 
-0.5.  The  <  rack  length  vs.  number  of  cycles  (a  vs  N)  data  for  the 
tests  are  shown  in  Figs.  3  through  5  (based  on  R  values). 

I 'sing  the  experimental  a  vs  X  data,  the  maximum  stress  intensity 
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Material  Characterization  Data.  Crack  Length  vs.  Life  to  Final  Fracture 
Maximum  Stress  -  10  ksi.  Load  Ratio  -  R=0.1 


Material  Characterization  Data.  Crack  Length  vs.  Life  to  Final  Fracture 
Maximum  Stress  -  10  ksi,  Load  Ratio  -  R=0.5 


Material  Characterization  Data.  Crack  Length  vs.  Life  to  Final  Fracture 
Maximum  Stress  -  10  ksi.  Load  Ratio  -  R=  -0.5 


factor  (K 


)  was  calculated  and  the  data  was  then  transformed  to 


max 

crack  maximum  growth  rate  vs.  maximum  stress  intensity  factor  (da/d.\ 
vs.  K  )  data  ns  inn  a  series  of  computer  programs  developed  by 
(Irimselv  [38),  which  are  based  on  the  ASTM-recommended  incremental 
polynomial-technique.  The  computer  programs  also  curve-fit  the  crack- 
growth  data  with  a  Walker  equation  [ 39 j .  The  (Irimselv  programs, 
called  the  Crack  Rate  Analysis  and  Walker  Equation  Solver  using  the 
Method  of  Least  Squares  (CRAWLS),  provide  a  major  portion  of  the  data 
and  equations  used  in  this  study  and  therefore  warrant  further 
d i scuss i on . 

The  CRAWLS  programs  transform  crack-growth  test  data  (a  vs  N) 
into  crack-growth  rates  (using  the  incremental  polynomial  technique) 
and  then  (its  those  growth  rate  data  with  a  Walker  equation. 

Conversion  of  crack-growth  data  (a  vs  X)  into  crack  growth  rate  data 
(da/dX  vs  K  max)  requires  differentiation  of  the  crack-growth  curve. 

I'll  i  s  differentiation  is  performed  using  a  method  which  fits  a  parabola 
to  sui  t  essi  ve  subsets  ot  seven  data  points  by  Lite  method  of  least 
squares  (see  e.g.  bv  liudar  el  al.  |4()|).  The  derivative  of  the  fitted 
curve  is  t  hc'n  calculated  for  the  middle  point  ot  each  subset .  This 
procedure,  combined  with  the  lames-Andersoii  backtracking  technique 
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Figure  h  :  Schematic  of  James-Anderson  Backtracking  Technique 


Note:  A  method  tor  solving  the  Walker  equation  for  C.,MtN  has  been 

proposed  by  Chang  et  al.  [41].  This  method  uses  a  least-squares  fit 
of  the  data  for  a  series  of  assumed  values  of  M.  Walker-equation 
plots  are  generated  and  chosen  based  on  which  best  fits  the  data. 

The  results  of  the  material-characterization  tests  are  given  in 
Figs.  7  through  9  for  the  three  values  of  the  R  ratio. 

Data  for  the  center-cracked-lension  specimens,  including  both 
positive  and  negative  stress  ratios,  were  input  into  Crimseiy's  CRAWLS 
programs.  The  positive-stress-ratio  data  can  be  represented  by  the 
following  best-tit  Walker  equation: 

_o  7  g  717Q7 

da/dN  =  1.11616  x  10  [K  f 1-R ) *  ]  ‘  (2) 

max 

where:  K  is  in  units  of  Ksi  per  in  and  da/dN  is  inch  per 

.  max 

c  yc  1  e 

The  best-fit  Walker  equations  for  the  negative  stress  ratios  are: 

_o  a  i  ]  | 

da/d.W  =  4.31344  x  10  (K  )  *  (3) 

max 

Four  fracture-toughness  tests  were  also  performed  for  the  2024- 
T331  aluminum.  These  tests  were  conducted  for  center-cracked  tension 
specimens  like  those  shown  in  Fig.  2.  The  fracture  toughness  was 
determined  by  a  combination  of  fatique  and  static  tensile  testing. 

Fhe  specimen  was  subjected  to  cyclic  loading  according  to  ASTM 
Standards  until  through-cracks  of  approximately  .003  inches  were 
initiated  on  both  ends  of  the  center  notch.  Then  the  specimen  was 
statically  loaded  bv  a  linearly  increasing  load  until  failure 


Figure  7  :  Material  Characterization  Data. 

Maximum  Stress  Intensity  Factor  vs.  Crack  Growth  Rate 
Maximum  Stress  -  10  ksi,  Load  Ratio  -  R=0.1 


Material  Characterization  Data. 

Maximum  Stress  Intensity  Factor  vs.  Crack  Growth 
Maximum  Stress  -  10  ksi,  Load  Ratio  -  R=0.5 


Material  Characterization  Data. 

Maximum  Stress  Intensity  Factor  vs.  Crack  Growth 
Maximum  Stress  -  10  ksi,  Load  Ratio  -  R=  -0.5 


occurred.  The  average  fracture  toughness  was  determined  using  die 
following  equation: 

K  =  J  -/  n.i  /sec(T7a/wX  (4) 

c  c 

Where:  0  c  »  stress  at.  failure 
a  =  crack  length 
w  =  specimen  width 

The  average  fracture  toughness  for  these  tests  was  found  to  be 
5  1.5  ksi  x  ./Tn  (48.68  MPa  x  /fir-). 

Standard  tensile  tests  using  a  "dogbone"  type  specimen  were 
performed  for  2024-T351  aluminum  to  determine  yield  strength.  The 
specimen  geometry  is  shown  in  figure  10.  The  average  yield  strength 
for  the  two  tests  performed  using  1/2-inch  thick  2024-1851  aluminum 
specimens  was  63  ksi  (913.59  MPa). 

2 . 2  Corner-Crack-at-a-llo  le  Tests 

Twelve  corner-crack-at-a-hole  tests  were  performed  using  2024- 
T351  aluminum  alloy  specimens.  The  typical  specimen  geometry  is  show 
in  Tig.  11.  The  specimens  were  precracked  at  the  same  constant 
amplitude  stress  levels  used  in  the  crack  growth  tests 
(  rJ  max  =  10  ksi).  Alter  prrcr,  uking,  the  maximum  stress  levels  used 
tor  testing  wore  IS  ksi  anu  1 5  ksi.  Two  posit ive  stress  ratios  ( R  = 

1 ) ,  1  and  0.3)  and  one  nog.it  ive  (K  =  -0.3)  stress  rat  io  were  considered 
i'iie  maximum  stress  I  eve*  I  ,  stress  rat  io,  and  initial  crack  size  and 
di.tpi  lot  ear  h  ol  the  I  we  I  vo  tests  are  given  in  Table  2.  Duplicate 


5625  inch 


tests  were  carried  out  for  each  of  the  six  load  conditions.  The 
results  were  averaged  to  eliminate  any  bias  that  might  arise  due  to 
any  possible  faulty  tests. 

To  evaluate  the  stress  intensity  factors  along  the  bore  of  the 
hole  for  the  2024-T351  aluminum  corner-crack-at-a-hole  tests,  the 
initial  crack  length  along  the  bore  of  the  hole  (a^)  must  be  known. 
Unlike  materials  such  as  PMMA,  aluminum  is  not  transparent  and 
currently  there  is  no  wav  to  measure  a^  during  the  test.  Therefore, 
to  obtain  values  for  ai  for  the  first  8  specimens,  a  marker  load  was 
applied  after  a  crack  was  observed  on  the  front  surface.  The  marker 
load  had  the  same  maximum  constant  amplitude  stress  levels  as  those 
used  in  the  tests,  but  the  minimum  stress  levels  were  increased  to  the 
stress  ratio,  R  =  0.85.  The  marker  load  produced  striations  on  the 
fracture  surfaces  which  could  then  easily  be  measured  with  an  optical 
microscope  after  the  completion  of  the  tests.  This  measurement 
yielded  the  value  of  ai  for  each  specimen.  During  testing  of  the  last 
four  specimens,  PC-9  through  PC-L2,  the  "Fax  Film"  method  was 
employed.  Fax  Film  is  a  thin,  mylar-tvpe  plastic  film  which  becomes 
very  soft  when  acetone  is  rubbed  on  it.  The  soft  film  is  forced  into 
the  small  initial  crack  and  then  removed,  a  crack  impression  is  left 
on  tile  film  and  with  the  aid  of  an  optical  microscope,  the  initial 
crack  length  can  be  measured.  This  method  is  very  accurate  if  the 
initial  crack  length  is  sufficiently  small;  it  also  requires  much  less 
time  and  effort  than  the  marker  band  technique. 
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The  corner-crack-at-a-hole  tests  tor  2024-T531  aluminum  specimens 
were  continued  until  failure  occurred.  Crack  length  measurements  were 
taken  throughout  the  test.  Crack  length  measurements  on  both  the 
front  and  back  surfaces  were  taken  at  random  intervals  during  the 
tests.  All  measurements  were  made  using  two  optical  microscopes 
mounted  on  the  fatique  machine,  one  at  the  front  surface  of  the 
specimen  anil  one  at  the  back  surface.  Using  the  crack  length 
measurements  and  the  number  of  fatique  cycles  between  each 
measurement,  the  crack  growth  rates  for  both  the  front  and  back 
surfaces  were  computed  using  methods  which  are  described  in  the  next 


chapter . 


3.0  EXPERIMENTAL  STRESS  INTENSITY  FACTORS 


3. 1 .Stress  Intensity  Factors  for  the  Front  and  Back  Surfaces 

Using  the  data  from  the  fatique  tests  experimental  stress 
intensity  factors  for  the  front  and  back  surfaces  were  obtained  using 
the  James-Anderson  backtracking  technique  [26].  The  technique  is 
shown  schematically  in  Fig.  6  and  may  be  briefly  outlined  (ref.  Opel 
[5],  p.23)  as  follows: 

1.  Crack-growth  curves  were  determined  at  the  back  surface  (c1 
vs  N)  and  the  front  surface  (c  vs  N)  using  raw  data  from  the 
corner-crack-at-a-hole  tests. 

2.  Using  the  crack  growth  curves,  crack-growth  rates  were 
determined  and  expressed  as  a  function  of  crack  lengths  (i.e., 
dc'/dN  vs  c'  and  dc/dN  vs  c).  Crack  growth  rates  were  determined 
using  the  standard  seven-point  polynomial  regression  technique 
discussed  earLier.  This  curve-fit  allows  for  the  calculation  of 
the  crack  growth  rate  at  any  crack  length.  Crack  length  vs. 
crack  growth  rate  are  plotted  for  the  front  surface  in  Figs.  12 
through  17  and  for  the  back  surface  in  Figs.  18  through  23. 

3.  The  stress  intensity  factors,  for  positive  stress  ratios, 
associated  with  given  crack  lengths  and  corresponding  crack 
growth  rates  were  found  using  the  Walker  equation.  The  Walker 
equation  form  for  the  positive  stress  ratios  was  found  to  be: 
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Figure  12:  2024-T351  Aluminum  Corner-Crack-at-a-Hole  Test  Data 

Front  Surface  Crack  Length  vs.  Crack  Growth  Rate. 
Maximum  Stress  -  18  ksi,  Load  Ratio  -  R=0.1 
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Figure  13:  2024-T351  Aluminum  Corner-Crack-at-a-Hole  Test  Data 

Front  Surface  Crack  Length  vs.  Crack  Growth  Rate. 
Maximum  Stress  -  18  ksi.  Load  Ratio  -  R=0.3 


2024-T351  Aluminum  Corner-Crack-at-a-Hole  Test  Data 
Front  Surface  Crack  Length  vs.  Crack  Growth  Rate. 
Maximum  Stress  -  15  ksi.  Load  Ratio  -  R=0.1 
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gure  16  :  2024-T351  Aluminum  Corner-Crack-at-a-Hole  Test  Data 

Front  Surface  Crack  Length  vs.  Crack  Growth  Rate. 
Maximum  Stress  -  13  ksi,  Load  Ratio  -  R=0.3 


Front  Surface  Crack  Length  vs.  Crack  Growth  Rate 
Maximum  Stress  -  15  ksi,  Load  Ratio  -  R=  -0.5 
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Figure  18  :  2024-T351  Aluminum  Corner-Crack-at-a-Hole  Test  Data 

Back  Surface  Crack  Length  vs.  Crack  Growth  Rate. 
Maximum  Stress  -  18  ksi,  Load  Ratio  -  R=0.1 
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Figure  l1):  2024-T351  Aluminum  Corner-Crack-at-a-Hole  Test  Data 

Back  Surface  Crack  Length  vs.  Crack  Growth  Rate. 
Maximum  Stress  -18  ksi,  Load  Ratio  -  R=0.3 
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igure  20 :  2024-T351  Aluminum  Corner-Crack-at-a-Hole  Test  Data 

Back  Surface  Crack  Length  vs.  Crack  Growth  Rate. 
Maximum  Stress  -  18  ksi,  Load  Ratio  -  R=  -0.5 


Figure  22  :  2024-T351  Aluminum  Comer-Crack-at-a-Hole  Test  Data 

Back  Surface  Crack  Length  vs.  Crack  Growth  Rate. 
Maximum  Stress  -  15  ksi.  Load  Ratio  -  R=0.3 


Figure  23  :  2024-T351  Aluminum  Corner-Crak-at-a-Hole  Test  Data 

Back  Surface  Crack  Length  vs.  Crack  Growth  Rate. 
Maximum  Stress  -  15  ksi,  Load  Ratio  -  R=  -0.5 


4.  0  ANALYTICAL  STRESS  INTENSITY  FACTORS 


Three  different  analytical  stress-intensity-factor  solutions  were 
considered  in  this  thesis.  The  first  solution  was  the  Newman-Raju 
three-dimensional  finite  element  solution  [9]  developed  for  corner 
cracks  at  a  hole.  The  second  solution  was  the  Grandt  linearization  of 
the  Bowie  equation  [10,  LI]  used  for  through  cracks  after  back 
surface  penetration  has  occurred.  Finally,  analytical  stress 
intensity  factors  were  calculated  using  Opel'  s  transition-region 
correction  factors  [5]  applied  to  the  Newman-Raju  solution  in  the 
corner  crack  region  and  to  the  Grandt/  Bowie  solution  for  a  through 
crack  after  back  surface  penetration.  After  the  stress-intensity- 
factors  were  calculated,  correlations  were  made  between  the 
experimental  results  and  the  Newman-Raju/  Grandt-Bowie  solutions 
without  transition  corrections  and  by  using  Opel'  s  transition  region 
correction  factors.  (These  results  will  be  used  in  a  later  section  to 
obtain  Life  predictions  in  order  to  assess  the  accuracy  of  the 
analytical  stress-intensity-factors. ) 

4 .  1  Newman-Raju  Solution  for  Corner-Cracks-At-A-Hole . 

The  Newman-Raju  solution  for  stress  intensity  factors  for  .corner 
cracks  at  holes  was  considered  because  of  its  accuracy  as  outlined  bv 
Heckel  and  Rudd  [36]  .  This  solution  was  also  chosen  because  it  is 
currently  one  ot  the  most  widely  used  solutions  in  the  aerospace 
industry.  The  equation  considered  was  based  on  an  initial  crack 


eccentricity  (a/c:  crack  depth  divided  by  crack  length,  see  Fig.l, 


greater  than  one).  The  equation  is: 

Kmax  ° max  Fch  (a/c’  a/t’  R/t’  c/b’  °)Fsh  (7) 

where : 

Q  =  1  +  1 . 464( c/a ) 1  • 6  5  (3) 

Fch  =  [M1+  M2(a/t)2  +  M3(a/tri  81  82  8 3  f(»  fw  (9) 

M,  =  /c7a  (1  +  0.04 (c/a))  (10) 

M?  =  0.2(c/aT  (11) 

>i3  =  -0.  1 1  (c/a)1*  (12) 

gl  =  l  +  [0.1+0.35(c/a)(a/t)2J(l-sin(t>)2  (13) 

g2  =  [ 1-0. 1 5 A  +  3.46A2  -  4.47A3  +  3. 52A4 ]/[ 1+0. 13A2  ]  (14) 

A  =  [  1  +  (c/r)cos(0.85<t>)  f 1  (15) 


g3  =  [  1 . 13-0.09(c/a )  J  [  1+0.  1  (  l-cos<j)  )2  ]  [0.8+0. 2(a/t)*25]  (16) 

£0  =  [(c/a)2  sin2(l)  +  cos2cJ)]*25  (17) 

f  =  [  sec(irR/2  b)sec (/a7t  (tt(  2R+c  ) )/( 4(  b-c  )+2c ) )  j • 5  (18) 

w 

F  h  =  [  ( 4/rr  +  ac/2tR)/(4/TT  +  ac/tR)]*5  (10) 


See  reference  9  for  details  of  this  solution. 


4.1.1  Boundarv-La ver  Effect  On  Stress  Intensity  Factors.  Hartranft 

and  Slh  [37]  suggested  that  the  stress  intensity  factors  in  a  very 

thin  "boundary- Layer"  near  the  intersection  of  the  crack-front  with  a 

free  surface  drop  off  rapidly  and  tend  to  zero  at  the  free  surface. 

This  idea  suggests  a  trend  opposite  to  that  of  Newman  and  Raju's  life 

prediction  solution  [9].  Their  solution  assumes  that  the  stress 

intensity  factors  tend  towards  infinity  at  the  surface.  Some 

reflection  on  the  problem  will  help  to  clarify  the  situation.  A  large 

stress  intensity  factor  would  lead  to  a  higher  value  of  da/dN,  which 

would  predict  a  shorter  life.  The  Newman-Raju  solution  has  been  known 

to  be  quite  conservative  [5,  6,  9]  i.e.  they  seem  to  over-estimate  the 

stress  intensity  factor  which  leads  to  a  prediction  of  shorter  lives 

than  those  observed  experimentally.  Newman  and  Raju  investigated 

Hartranft  and  Sih's  suggestion  using  finite  element  techniques  and 

found  that  the  stress  intensity  factors  tended  toward  zero  at  the  free 

surface,  as  proposed.  This  implies,  theoretically,  that  no  crack 

would  grow  at  the  surface  of  the  specimen  (i.e.,  if  K  approaches 

max 

zero,  da/dN  approaches  zero,  da/dN  approaches  zero).  Newman  and  Raju 
subsequently  investigated  the  effect  of  the  boundary  layer  by 
evaluating  stress  intensity  factors  along  an  imaginary  surface  at  an 
angle  d>  from  the  free  surfaces  (see  Fig.  1).  A  graph  showing  the 
trend  of  the  Newman-Raju  correction  and  the  result  of  evaluating  the 
stress  intensity  factors  at  an  angle  d>  from  the  free  surface  is  shown 
in  Fig.  24.  It  can  easily  be  seen  that  at  the  free  surface  (i.e., 


where  2  $  /  tt  =  l),  there  is  some  inconsistency.  If  the  angle  is 
increased,  the  solutions  become  more  consistent.  (Mewman  and  Raju  did 
find  that  beyond  a  value  of  <j>  >  .4tt,  there  was  no  difference  between 

the  two  solutions,  see  Fig.  24). 

Recently,  an  ASTM  Task  Group  was  established  for  the  purpose  of 
improving  current  methods  for  calculating  stress  intensity  factors. 

The  task  group  has  been  sponsoring  an  annual  stress  intensity  factor 
solution  "round-robin".  The  participants  provide  stress  intensity 
factor  solutions  and  these  are  evaluated  by  comparing  them  to  other 
participant  solutions  for  the  purpose  of  finding  better  solutions. 

Some  participants  have  used  this  idea  of  evaluating  the  stress 
intensity  factor  along  a  surface  at  an  angle  of  $  (  $  ,  is  the 
parametric  angle  and  is  equal  to  0°  at  the  front  surface  and  90°  at 
the  base  of  the  hoLe),  most  notably  Chung  [43],  and  have  met  with  a 
great  deal  of  success.  An  appropriate  value  for  the  angle  0,  based 
on  results  from  the  round-robin  competit'1  on ,  is  <J>  =  10°. 

Evaluation  of  the  stress  intensity  factors  along  a  surface  at 
<+)  =  10°  (see  Fig.  1  for  geometry),  has  lead  to  improvement  in  accuracy 
of  life  prediction  using  the  \ewman-Ra ju  corner-crac k-af -a-hole 
solution.  Life  Predictions  until  back  surface  penetration  (the  only 
region  where  rf>  has  an  effect)  improved  an  average'  of  20%  for  2H24- 
T331  aluminum  and  30%  for  7073-1031  aluminum  using  the  Ctu  orrec ted 
N'ewman-Ro  ju  sollution  (see  Tables  3  and  4,  respectively).  The 
improvement  in  life  prediction  is  quite  significant  and  therefore,  the 
solutions  for  2024-1531  and  7073-1031  aluminum,  studies  in  this  thesis 
will  include  the  angle  =  10°  in  all  calculations. 


re  on 


90*  7.b) 


4 .  2  Grandt  s  Linearization  of  Bowie  s  Solution  for  Through  Cracks. 
The  Grandt  linearization  of  the  Bowie  equation  described  by  Engle  [11] 
was  considered  for  through  cracks  from  back  surface  penetration  until 
final  fracture.  The  equation  considered  was: 


K  =0  /ttc  F  (c/r)f 
max  max  gr  w 


where : 


Details  for 


F  =  .6762062  =  [ . 873301 5/( . 3245442(c/R )  ] 

gr 

f  is  given  by  equation  (18). 
this  solution  are  given  in  reference  11. 


(20) 


4.3  Opel's  Correction  Factors.  In  the  third  of  the  solutions 

considered,  in  his  thesis  Opel  introduced  transition  region  correction 

factors  for  the  Newman-Raju  and  Grandt/Bowie  solutions.  Opel 

separated  the  transition  region  into  two  different  regions.  The 

beginning  of  the  first  transition  region  was  where  a/t  =  0.75  and 

ended  where  a/t  =  1.0  (see  Fig.  25  for  details).  From  the  initial 

crack  length  to  where  a/t  =  0.75,  the  Newman-Raju  solution  was  used 

without  correction  factors.  From  a/t  =  0.75  to  a/t  =  1.0  Opel  applied 

the  correction  factor  F  to  the  Newman-Raju  solution.  The  new 

tr 

stress-intensity-factors  (K)  new  were  found  using  an  equation  of  the 
form : 

(K)  .  „ 


=  (K)m  o  X  (F) 


(21  ) 


a/t^O.75 


F  ,  OpeL's  correction  factor,  was  found  using  the  James-Anderson 
backtracking  technique  and  a  polynomial  regression  curve  fit,  and  was 
determined  to  be  of  the  form 

(Ftr}b  =  ^I-573924  +  l .468068(a/t)-l .275787(a/t)2]  , 
for  the  base  of  the  hole  and 

(F  )  =  1/ [ 1.31687  -  . 362863(a/ 1 )  ] ,  a<t 

for  the  front  surface. 

Transition  correction  factors  were  also  developed  for  the  region 
from  the  back  surface  penetration  (a/t  =  1.0)  to  the  end  of  the 
transition  region  where  c/r  =  2.5.  The  correction  factor  for  the  back 
surface  was  found  to  be: 

(F  )b  =  1/ [ . 191601  +  . 724961 ( c ' /R )  -. 1 58451 ( c ' /R )2 J 

The  equation  used  for  the  stress-intensity-factors  on  the  back  surface 
transition  region  was: 

(K).  =  /cHT Fp  ( c  '  / R )  Fw(  F  ). 

bs  B-0  tr  bs 

where  FD  .,(c'/r)  is  the  C.randt  curve  fit  (10J  to  Bowie's  tabular  hole 

iJ  vi 

.orrection  factor  [12].  fw  is  the  iinite  width  correction  used  by 
Newman-Raju  [ ')  ] .  The  Grand t  linearization  of  the  Bowie  solution  was 
used  for  stress  intensity  factors  for  the  region  from  the  end  of  the 
transition  region  to  final  fracture. 


5.0.  DEVELOPMENT  OF  TRANSITION  CRITERIA 


As  the  crack  propagates  towards  back  surface  penetration,  it  retains  a 
quarter-elliptical  shape  with  varying  aspect  ratio  until  a  =  t.  (See 
Fig.  1.)  Therefore,  no  correction  factor  should  be  necessary.  This 
idea  contradicts  the  work  by  Opel,  but  there  are  good  reasons  that  can 
be  justified  analytically. 

Opel  developed  transition  correction  factors  in  this  region 
primarily  to  improve  on  Newman  and  Raju's  conservative  life  prediction 
solution  for  corner  cracks  at  a  hole.  The  discussion  of  the  boundary 
layer  effects  in  the  previous  section  suggested  by  Hartranft  and  Sih 
[37]  and  tested  by  Newman  and  Raju  [9]  provide  a  means  by  which  to 
improve  the  conservative  Newman-Raju  corner  crack  solution. 

The  Newman-Raju  solution  is  a  function  of  the  angle  ,  among 
other  parameters.  (See  Fig.  1.)  Opel's  correction  factors  tor  stress 
intensity  and  life  predictions  are  based  on  <t>  =  0 0  and  90°, 
corresponding  to  the  front  surface  and  the  surface  along  the  base  of 
the  hole  respectively.  Hartranft  and  Sih  suggest  that  calculation  of 
stress  intensity  factors  along  an  imaginary  surface  at  the  interior  of 
the  specimen  would  lead  to  possibly  more  accurate  life  predictions. 
This  idea  seems  to  be  valid  based  on  finite  element  solutions  made  bv 
Newman  and  Raju.  In  fact,  life  predictions  are  improved  drastically 
i  f  <J)  is  10°  and  80°  which  correspond  to  imaginary  surfaces  near  the 
front  surface  and  near  the  surface  at  the  base  of  the  hole 
respectively.  Tables  5  and  6  show  comparisons  of  solutions  for  2024- 
T351  and  7075-T651  aluminums  respectively.  The  first  solution  is  the 


Using  Newman-Ra  ju  ( Uncorrer ted )  fur  7075-T6‘il 


uncorrected  N'ewman-Raju  solution  with  0  =  0°  and  90°  and  the  second  is 
the  uncorrected  Newman-Raju  solution  with  <1  =  10°  and  80°  .  Comparing 
these  two  solutions  for  Life  untiL  back  surface  penetration  (a=t),  it 
can  easily  be  seen  that  using  0  =  10°  and  80°  improves  the  Newman- 
Bowie  model  enough  to  warrant  the  exclusion  of  all  correction  factors 
used  for  life  until  back  surface  penetration. 

5.1.2.  Back  Surface  Penetration  Until  Final  Fracture.  Opel  chose  to 
use  transition  region  correction  factors  in  this  region  based  on  front 
and  back  surface  crack  lengths,  normalized  by  the  hole  radius  (R),  :/R 
and  c'/R  respectively.  These  transition  region  correction  factors 
were  applied  from  back  surface  penetration  (a  =  t)  to  a  point  where 
c/R,  c'/R  =  2.5.  (See  Opel  for  discussion  on  development  of  the  end 
point  for  the  transition  region.)  The  Opel  transition  correction 
factors  produced  favorable  results  towards  increasing  the  accuracy  of 
the  Orandt  linearization  of  the  Bowie  solution  for  through  the 
thickness  cracks.  However,  when  examined  closely,  it  turns  out  that 
these  corrections  are  highly  dependent  on  specimen  thickness.  This 
dependence  can  easily  be  seen  by  considering  an  example  (see  Fig.  26). 
The  specimens  considered  by  Opel  were  fabricated  from  7075-T651 
aluminum  with  thickness  t  =  .25  inches  and  a  hole  radius  R  =  .125 
inches.  After  the  crack  propagates  through  the  back  surface  at  a  =  t, 
Opel  applied  the  transition  correction  factors,  based  on  c/R  and  c'/R, 
to  the  Grandt-Bowie  solution  with  favorable  results.  Opel  applied 
correction  factors,  based  on  specimens  with  t  =  .25,  to  the  region  a  = 


(or 


c'/R  =  O'  to  c'/R  =  2.5. 


At  c'/R  =  2.5,  the  crack  was 


considered  a  through-the-thickness  crack  and  after  that  point  no 
corrections  were  necessary  to  modify  the  Grandt-Bowie  solution. 

(Note:  c'/R  =  2.5  corresponds  to  c'  =  .3125  inches  and  at  this  point 
c  =  .3125  inches  as  well.)  If  one  were  to  consider  a  thicker 
specimen,  say  t  =  1.0  inch,  with  the  same  hole  radius  (R  =  .125"),  the 
crack  would  not  have  reached  back  surface  penetration  (a  =  t)  when  c/R 
=  2.5  (i.e.  c'  =  0).  Therefore,  applying  surface  correction  factors 
at  c/R  =2.5  for  the  thicker  specimen  would  be  invalid.  Clearly, 
correction  factors  based  on  t  =  .25  and  dependent  on  c/R,  c'/R  are 
only  valid  for  specimens  .25"  thick  and  new  correction  factors  must  be 
developed  for  different  thicknesses.  Developing  new  correction 
factors  for  each  thickness  would  provide  favorable  results  but  would 
be  much  too  tedious  a  task,  involving  experiments  and  analytical 
corrections  each  time  a  life  prediction  is  necessary.  Therefore,  the 
development  of  correction  factors  that  are  independent  of  specimen 
thickness  is  highly  desirable. 

Rudd  [ 42 Jsuggested  that  the  transition  region  after  back  surface 
penetration  be  corrected  by  transition  correction  factors  based  on  the 
shape  of  the  crack  as  it  propagates  from  a  =  t  to  a  complete  through- 
the-thickness  crack.  The  shape  of  the  crack  can  be  dealt  with  by 
understanding  how  the  front  surface  crack  length  (c)  and  the  back 
surface  crack  length  (c')  vary  with  each  other.  When  c'  =  c  (or  c'/c 


=1),  the  crack  lengths  on  both  surfaces  are  equal  and  a  true  through- 


the-thickness  crack  exists  (See  Fig.  1.)  If  transition  correction 
factors  are  based  on  c'/c,  then  the  correction  factors  become 
th  Lckness-independent . 

5.2.  Anal ytical/ Experimental  Correlations 

Correlations  were  made  of  the  experimental  stress  intensity 
factor  ranges,  using  a  form  of  the  Walker  equation,  with  the 
analytical  predictions,  based  on  the  Grandt-Bowie  solution,  for  2024- 
T351  and  7075-T65!  aluminum  for  the  region  from  back  surface 
penetration  until  the  final  fracture  region.  Experimental  stress 
intensity  factors  were  obtained  using  the  James-Anderson  backtracking 
technique,  based  on  duplicate  specimens  for  each  of  6  different  load 
conditions.  The  procedure  used  to  obtain  the  appropriate  correlations 
is  outlined  next: 

1)  Using  the  James-Anderson  backtracking  techniques, 
experimental  data  in  the  form  of  c  vs  N  (and  c'  vs  M  for  back 
surface)  were  converted  to  c  vs  dc/dN  (and  c'  vs  dc'/dN). 

2)  Results  from  1)  above  for  each  duplicate  set  of  specimens 
were  curve-fitted  using  an  ASTM-recommended  polynomial 
regression  scheme.  The  data  were  "best  fitted"  to  5th-order 
equations.  This  provided  dc/dN  and  dc'/dN  equations  as 
functions  of  c  and  c'  for  each  set  of  duplicate  specimens. 


3)  The  experimental  stress  intensity  factor  (I(„)  was  calculated 

h. 

from  back  surface  penetration  to  final  fracture  using  the 
Walker  equation. 

4)  The  analytical  stress  intensity  factor  (K^)  was  calculated 
from  back  surface  penetration  to  final  fracture  using  the 
Grandt-Bowie  solution. 


5)  Results  are  plotted  in  the  form  of  c'/c  vs  K^/!(g  (See  Figs. 
27  to  30.) 

The  correlations  show  that  at  c'/c  =  1.0,  K^/K^,  approaches  a  constant 
Therefore,  the  end  of  the  transition  region  is  where  c'/c  =  0.9  and 
the  crack  can  be  approximated  as  a  through-crack.  This  result  can 
favorably  be  compared  to  Johnson  [34]  who  defined  the  end  of  the 
transition  region  to  be  c'/c  =  .9  for  life  prediction  approximations. 
5.3.  Correction  Factor  Development 

5.3.1.  Back  Surface  Penetration  Until  Final  Fracture.  As  noted  in 
the  previous  section,  correlations  tend  to  become  constant  at  K^/K^ 
when  c'/c  =  0.9.  Fig. 31  shows  the  correlations  for  the  back  surface 
data.  The  polynomial  regression  curve  fit  for  this  data  was  found  to 
be 

( Ft  r ) ,  =  [ l/.609633-.689614(c’/c)+l .23275(c’/c)2  ] 

bs 

and  it  is  plotted  in  Fig. 31 .  It  was  found  that  a  quadratic  curve  fit 


of  the  2024-T351  Aluminum  Test  in  the  Region  of  Interest 


Correlation  of  Stress  Intensity  Factors  Along  the  Front 
Surface  of  the  7075-T651  Aluminum  Test  in  the  Region  of  Interest 


of  the  2024-T351  Aluminum  Test  in  the  Region  of  Interest 


of  the  7075-T651  Aluminum  Test  in  the  Region  of  Interest 


(Correction  Factor  Plotted  in  the  Transition  Region. 


i  ••--sil r  i:,  t  g* "  ;ii  correlation  for  the  data.  Correlations  did  not 
>vo  si  ; : i  i !  i i  intly  when  higher  order  curve  fits  were  attempted, 
i'll  is  correct  ion  is  applied  to  the  Grandt-Bowie  solution  after 
buck  surface  penetration  until  final  fracture.  This  correction,  when 
applied  to  the  correlation  data,  essentially  brings  all  the  data 
points  to  a  constant  value  of  K^/Kg  =  1.0.  The  effect  would  be  to 
"speed  up"  the  growth  of  the  crack  on  the  back  surface,  and  "slow 
down"  the  growth  of  the  crack  on  the  front  surface.  Therefore  more 
accurate  life  predictions  will  result  when  the  correction  factors  are 
used  to  modify  the  stress  intensity  factors  on  the  front  and  back 
surfaces . 

"j.3.2.  Front  Surface . 

In  a  similar  manner,  the  front  surface  correction  factor  based  on 
the  polynomial  regression  curve  tit  is: 

(Ftr}fs  =  t1/]  *24"')7  -  •  1  '/O  I 

and  it  is  plotted  in  Fig. 32.  Again,  it  was  found  that  a  linear  curve 

fit  resulted  in  a  good  correlation  and  did  not  improve  significantly 
when  higher  order  curve-fits  were  attempted. 

The  correction  factors  for  the  front  and  back  surfaces  were 
substituted  into  the  Grandt-Bowie  equation  a.,  follows: 

^max^cor  ~  ^max  V,-B^Ftr  ^bs 
^'max^cor  -  ^ir.ax ^G-B'  ^ F  t  r  ^  f  s 

where  (K  )  is  the  corrected  stress  iniensitv  actor,  (K  ),.  p  is 
max  cor  max  G-b 

the  Grandt-Bowie  solution  and  (Ftr)  ,  ( Ft r )  are  the  hack  ini  front 

1)  s  I  s 

surface  correction  factors  respectively. 


Figure  32:  Superposition  of  Correlations  of  Stress  Intensity  Factors  Along 
the  Front  Surface  for  2024-T351  and  7075-T651  Aluminum  Tests. 
(Correction  Factor  Plotted  in  the  Transition  Region.) 


In  summary,  the  growth  of  a  crack  as  it  grows  from  a  corner- 
crack-at-a-hole  to  final  fracture  can  be  separated  up  into  three 


phases : 

1)  Corner  crack  at  a  hole  until  back  surface  penetration  (a  = 
t);  this  phase  is  dealt  with  by  evaluation  for  stress 
intensity  factors  along  "imaginary"  surfaces  near  the  top 
surface  and  the  surface  at  the  base  of  the  hole  to  preclude 
effects  of  fabrication  (cold-rolling,  hole  drilling,  etc.). 
No  other  corrections  are  necessary  to  the  Newman-Raju 
solution . 

2)  Back  surface  penetration  (a  =  t)  until  through-the-thickness 

crack  (c'/c  =  0.9);  this  is  dealt  with  by  applying  front-and 

back-surface  correction  factors,  (F  )c  and  (F_  ), 

tr  fs  tr  bs 

respectively,  to  the  Grandt-Bowie  solution  for  through-the- 
thickness  cracks. 


3)  Through  the  thickness  crack  until  final  fracture;  use  only 
the  Grandt-Bowie  solution,  no  corrections. 


6.0  Evaluation  of  Transition  Criteria. 


6 . 1  Life  Prediction  Models  Used 

An  evaluation  is  made  of  the  accuracy  of  the  transition  criteria 
developed  in  this  study  as  well  as  that  of  an  instantaneous  criterion 
[11],  Opel's  criterion  [5],  a  criterion  developed  by  Brussat  et  al. 

[8]  and  one  developed  by  Collipriest  and  Ehret  [7].  Life  predictions 
are  made  using  these  criteria.  These  predictions  are  then  compared 
with  each  other,  as  well  as  with  the  corner  crack  at  a  hole  test 
results  for  7075-T651  and  2024-T351  aluminum  previously  described. 

Life  predictions  were  made  for  corner  crack  at  a  hole  specimens 
subject  to  constant  amplitude  loading.  These  life  predictions  were 
based  on  material  characterization  data  obtained  for  2024-T351  and 
7075-T651  aluminum.  This  characterization  data  was  used  to  derive  the 
proper  form  of  the  Walker  equation  for  each  type  of  aluminum,  see 
equations  2  through  3  for  2024-T351  aluminum  and  he  following 
equations  were  used  for  7075-T51  aluminum  (from  Opel  [5]): 

da/dN  =  3.2624  x  10~9[ ( i-R ) ’ 5K  ]3,3908  (23) 

max 

for  positive  stress  ratios  and  for  negative  stress  ratios: 

da/dN  =  1.29  x  10_8[K  ]2,89  (2A) 

max 

Using  the  Walker  equations,  the  crack  growth  rate  was  calculated  and 
then  the  new  crack  length  was  found  by  integrating  after  each  cycle. 
Subsequent  crack  lengths  were  found  by  adding  the  change  in  crack 
length  to  the  previous  crack  length. 
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Three  different  lives  were  considered  for  each  life  prediction 
model:  (a)  from  the  initial  crack  sizes  specified  in  Tables  2  and  7 

to  back  surface  penetration  (a=t),  (b)  from  back  surface  penetration 


to  final  fracture  and  (a=t -►  0=0^ ) ,  (c)  total  life  (C  =  C  ).  The 
instantaneous  transition  criterion  predictions  were  made  using  the 
Newman-Raju  corner  crack  at  a  hole  solution  [9]  until  back  surface 
penetration  occurred  (i.e.  a  =  t).  The  Grandt  linearization  of  the 
Bowie  solution  for  through  the  thickness  cracks  [10]  was  then  used, 
assuming  an  initial  back  surface  crack  length  equal  to  the  front 
surface  crack  length  at  the  time  of  back  surface  penetration.  In 
addition  to  the  transition  criteria  developed  in  this  study,  the 
Brussat,  Opel  and  Coll ipriest-Ehret  criteria  were  applied  according  to 
procedures  outline  in  [8,  5,  and  7]  respectively. 

6.2  Results. 

6.2.1.  Life  Predictions  Based  on  Corrected  Stress  Intensity 
Factors.  (Note:  Experimental  results  were  obtained  bv  Heckel  &  Rudd 
[25]  and  analytical  predictions  made  by  Opel  [5]  for  a  7075-T651 
aluminum.  Experimental  lives  for  2024-T351  and  7075-T651  aluminum  are 
shown  in  Tables  8  and  9,  respectively.) 

Analytical  life  predictions  were  made  for  the  six  2024-T351 
aluminum  test  conditions  presented  in  Table  2  and  the  seven  7075-T651 
aluminum  test  conditions  presented  in  Table  7  using  the  criteria 
developed  in  this  study.  Each  experimental  life  shown  in  these  tables 
represents  the  average  of  duplicate  specimens.  The  life  predictions 


based  on  the  corrections  developed  in  this  thesis  are  shown  in  Tables 
10  and  11  for  2024-T351  and  7075-T651  aluminum,  respectively. 


6.2.2  Life  Predictions  Based  on  Other  Models. 

6.2.2. 1  Newman-Bowie .  The  Newman-Bowie  (instantaneous) 
method  was  used  as  presented  by  Engle  [11]  to  predict  the  fatique  life 
for  the  specimens  referenced  in  the  previous  section.  The  important 
consideration  for  this  model  is  the  fact  that  it  ignores  the 
transition  region  altogether.  It  assumes  a  through-the-thickness 
crack  once  back  surface  penetration  has  occurred.  The  results  are 
presented  in  Table  12  for  2024-T351  and  Table  13  for  7075-T651 
aluminums . 

6. 2. 2. 2  Opel .  The  development  of  Opel's  correction  factors 
are  discussed  in  section  4.3.  (See  Opel  [5]  for  more  details).  His 
transition  region  correction  factors  for  the  base  of  the  hole  were 
derived  from  PMMA  data  and  based  on  $  =  0°  and  90"  .  The  transition 
correction  factors  for  the  front  and  back  surfaces  are  based  on 
analytical  and  experimental  stress  intensity  factors  as  a  fu;nction  of 
c/R  and  c'/R  for  7075-T651  aluminum.  Life  predictions  using  these 
transition  region  correction  factors  are  shown  for  2024-T351  aluminum 
and  7075-T651  aluminum  in  Tables  1  and  14,  respectively. 

6. 2. 2. 3  Brussat .  The  Brussat  method  [8]  was  developed  for 
quarter-circular  cracks  at  holes,  therefore  the  initial  crack  length 
used  was  the  average  of  the  initial  crack  lengths  along  the  base  of 
the  hole  and  along  the  front  surface.  The  procedure  used  for  the 
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Brussat  Lite  predictions  are  the  same  as  outlined  in  section  A,  except 
for  calculating  the  stress  intensity  factors  where  Brussat  used  the 


following  equation: 


K  =0  Ac  f  f 
max  max  w  tr 


where  f  =  1  -  (0.2886/(1  +  2(c/t)2- 


and  f  is  the  finite  width  correction  by  equation  (18). 


i’he  results  of  the  Brussat  predictions  are  shown  in  Table  15  for  2024- 
!B5i  aluminum  and  Table  16  for  7075-T651  aluminum. 

6. 2. 2. 4  Coll ipr iest-Ehret .  Collipriest  and  Ehret  [7] 
suggested  transition  criteria  based  on  the  front  and  back  surface- 
crack  lengths.  Similar  to  the  method  developed  in  this  study, 
Collipriest  and  Ehret  did  not  use  any  correction  until  after  back 
surface  penetration.  The  Newman-Bowie  prediction  method  was  used  as  a 
basis  for  this  life  prediction.  The  Newman-Raju  corner  crack  at  a 
hole  solution  was  used  to  calculate  stress  intensity  factors  until 
back  surface  penetration  (a  =  t).  Then  an  imaginary  crack  length 
along  the  base  of  the  hole  was  used  of  the  form: 


a’  =  T(  1  -  ((’  /Cr)2]  *5 


where  Cy  =  front  surface  crack  length 

C  =  back  surface  crack  length 
\) 

T  =  specimen  thickness 


0  max 


This  crack  length  was  used  in  the  Newman-Raju  solution  for  the  stress 
intensity  factors  along  the  front  surface. 

The  stress  intensity  factors  along  the  back  surface  were  calculated 

using  the  Grandt-3owie  through  crack  solution  of  the  form: 

Kr  =0  max  /  tt  c  F  (c/R)  F  f 
Lmax  gr  w  cec 

f  =  [  l/( 1-(1-(C  /C  )2))]'5 
cec  b  t 

and  F  (c/R)  and  f  are  given  by  equations  22  and  18. 
gr  w  °  n 

When  the  back  surface  crack  length  equals  the  front  surface  crack 
length,  only  the  uncorrected  Grandt-Bowie  equation  was  used  until 
final  fracture.  (See  Fig. 33  for  details  of  imaginary  crack  length.) 
The  results  of  the  Coll ipriest-Ehret  models  are  shown  in  Tables  17  and 
18  for  2024-T351  and  7075-T651  aluminum,  respectively. 

6.3.  Discussion  of  Results 

The  model  developed  in  this  thesis  is  the  most  accurate  and 
versatile.  The  average  total  lives  for  2024-T351  and  7075-T651 
aluminum  are  shown  in  Tables  19  and  20.  All  predictions,  except  for 
Opel's  and  Brussat's,  have  been  made  using  the  parametric  angle 
(  S>  =  10°  ,  80°  )  corrections.  The  new  model  shows  an  average 
improvement  of  4%  over  the  instantaneous  model.  This  improvement  is 
significant,  considering  that  the  transition  region,  where  the  new 
model  differs  from  the  instantaneous  model,  represents  only  a  small 
portion  of  the  specimen's  total  life.  The-  new  model  also  shows  an 
improvement  of  3%  over  the  Opel  model.  However,  it  should  also  be 


[jointed  out  that  the  Opel  model  is  limited  since  the  transition  region 
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correction  factors  it  uses  were  derived  specifically  for  specimens 
0.25  inches  thick,  and  must  be  rederived  for  other  thicknesses. 

The  Brussat  model  seems  to  be  more  accurate  than  the  new  model , 
but  a  closer  look  will  reveal  that  this  is  not  true.  As  shown  in 
fable  16,  the  Brussat  model  is  unconservative  for  three  out  of  seven 
tests.  However,  an  examination  of  Table  16  reveals  that  the  Brussat 
predictions  for  life  until  back  surface  penetration  very 
unconservative,  while  those  for  the  region  from  back  surface 
penetration  until  final  fracture  are  conservative.  This  accounts  for 
the  good  average  total  life  predictions.  Table  15  shows  that  the 
Brussat  predictions  for  2024-T351  aluminum  which  is  more  ductile  than 
the  7075-T651  aluminum  are  much  more  accurate.  Therefore,  these 
results  seem  to  indicate  that  caution  should  be  exercised  in  the 
application  of  the  Brussat  model  to  materials  which  are  relatively 
brittle.  On  the  other  hand,  the  model  put  forth  in  this  investigation 
seems  to  yield  conservative  predictions  for  both  relatively  ductile 
and  brittle  aluminum  alloys. 

The  Col  1 i pr i est-Ehret  model  was  found  to  be  limited  when  applied 
to  thin  specimens.  As  previously  discussed,  this  model  calculates  the 
stress  intensity  factors  for  the  front  surface  after  back  surface 
penetration  using  an  imaginary  crack  length  in  the  Newman-Ra ju  corner- 
crack  solution.  Due  to  the  form  ot  the  equation  used  to  calculate  t he 
imaginary  crack  length,  a  relatively  large  initial  back  surface*  crack 
length  (of  the  order  of  0.  I  inch)  must  be*  assumed.  When  this  is  done 
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fur  chin  specimens,  a  good  deal  of  the  transition  region  is  ignored, 
and  therefore  the  predictions  become  almost  identical  to  those  ot  the 
instantaneous  model.  Therefore,  there  seems  Lo  be  no  advantage  in 
applying  this  model  to  relatively  thin  specimens  instead  of  the 
instantaneous  model. 
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This  investigation  showed  that  although  Opel's  model 
yielded  more  accurate  life  predictions  for  a  new,  more  duc¬ 
tile  material  (2024-T351  aluminum),  it  was  very 
conservative  in  predicting  total  life.  The  new  model 
proposed  in  this  work  was  shown  to  be  considerably  more 
effective. 

The  new  model  was  developed  by  first  considering  the 
Hartranft  and  Sih  suggestion  to  modify  the  Newman-Raju 

solution  with  a  parametric  angle  (j>  =  10“  and  80°  insteadof 
0°  and  90*  traditionally  used.  This  correction  eliminated 
the  surface  boundary  effects  and  resulted  in  significant 
improvements  in  the  life  prediction  obtained  from  all 
models  that  were  investigated.  As  a  result  of  using  f>  =  10 
and  80°,  no  transition  correction  factors  were  necessary 
for  the  bore  of  the  hole  in  the  life  predictions  from 
corner-crack  to  back  surface  penetration. 

The  Newman-Raju  corner-crack  and  the  Grandt 
linearization  of  the  Bowie  through-crack  solutions  were 
found  to  be  accurate,  but  conservative  in  predicting  the 
stress  intensity  factors  for  a  corner-crack  emanating  from 
a  hole  as  it  grows  until  final  fracture.  However,  neither 
of  these  solutions  considers  the  transition  region  where 
the  crack  propagates  neither  as  a  corner-crack  or  a 
through-crack.  This  transition  region  was  determined  to 


begin  at  back  surface  penetration  and  end  when  the  back 
surface  crack  length  equals  the  front  surface  crack  length 
(i.e.,  c'  =  c).  At  this  point  the  crack  becomes  a  true 
through- the- thickness-crack.  Transition  correction  factors 
based  on  correlations  of  analytical  and  experimental  stress 
intensity  factors  were  derived  as  a  function  of  c'/c. 
Using  transition  factors  based  on  c'/c  instead  of  c/R,  as 
Opel  did,  eliminated  the  th i ckn e s s- d e p e nd e n ce  of  the 
transition  factors  and  as  a  result  they  are  more  versatile. 

The  1 if e- pr e di c t i on  model  developed  in  this  thesis 
proved  to  be  the  most  accurate  and  versatile.  It  accounts 
for  the  transition  region  while  the  Instantaneous  model 
does  not.  It  uses  th i c k ne s s- i n de pe n d e n t  transitions 
correction  factors  while  the  Opel  transition  factors  have 
to  be  re-derived  when  the  thickness  changes.  The  new  model 
yields  consistently  conservative  predictions,  while  the 
Brussat  model  is  very  un co n s e rv a t iv e  for  predicting  life 
until  back  surface  penetration  and  is  conservative  for  life 
predictions  from  back  surface  penetration  until  final  frac¬ 
ture.  Finally  it  uses  correction  factors  which  apply  over 
any  thickness  while  the  Co  1 1 i p r i e s t- Eh  re t  mode]  seems  to 
assume  a  relatively  large  initial  back  surface  crack  length 
which  makes  it  inappropriate  for  thin  specimens. 

This  study  has  demonstrated  that  the  proposed 
corrections  to  the  currently  used  stress  intensity  factor 
solutions  have  resulted  in  better  life  predictions  for 


constant  amplitude  loading.  The  life  prediction  model 
developed  in  this  study  showed  a  three  or  four  percent 
improvement  over  the  uncorrected  life  prediction  model. 
Even  though  this  does  not  seem  like  a  sizable  improvement, 
one  must  keep  in  mind  the  specimen  thickness  considered  was 
only  0.25  inches.  Therefore,  the  transition  region  is  a 
very  small  part  of  the  total  life  of  the  specimen.  Further 
investigation  using  this  prediction  model  and  thicker 
specimens  should  be  undertaken.  Thicker  specimens  would 
have  a  larger  transition  region  and  the  new  model  should 
result  in  greater  improvements  in  the  accuracy  of  the 
predictions.  Applications  of  this  model  to  lugs,  for 
instance,  may  prove  useful.  Lugs  have  geometries  such  that 
there  are  thick  specimen  depths  and  narrow  specimen  widths. 
Therefore  the  transition  region  would  account  for  a 
significant  part  of  the  total  life  after  back  surface 
penetration.  Other  follow-on  topics  should  include  the 
investigation  of  the  effect  of  different  parametric  angles, 
(j>  .other  than  j  =10*  and  80*  as  well  as  applying  the  new 
model  to  different  types  of  material,  such  as  steel. 
Should  the  correction  factors  developed  in  this  thesis 
prove  to  be  as  effective  as  they  seem  to  be,  they  should  be 
considered  for  incorporation  into  the  Air  Force  Damage 
Tolerance  Design  Handbook  [44]. 


Computer  Programs  Used 
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There  were  six  major  computer  programs,  written  in 

Fortran  77,  used  for  this  study.  A  short  discussion  of 

what  each  program  does,  followed  by  complete  program 

listings  are  included  in  this  appendix.  The  input  variables 

for  the  life  prediction  programs  are: 

A  -  Initial  crack  length 

C  -  Initial  crack  depth 

R  -  Radius  of  hole 

S  -  Maximum  stress 

T  -  Thickness  of  specimen 

W  -  Width  of  specimen 

ASTOP  -  Non- f un cti onal  parameter 

RS  -  Stress  ratio 

Cl, N, FAC  -  Walker  constants  (mat'l  dependant) 

A.  Newman-Bowie  ( Un corre cted ) 

Also  known  as  the  instantaneous  model,  this  is  the 
program  upon  which  all  of  the  other  life  prediction 
programs  are  based  (except  Brussat's  model).  It  uses  a 
Newman-Raju  corner-crack  solution  for  stress  intensity 
factors  until  back  surface  penetration.  After  back  surface 
penetration,  the  program  uses  a  Grandt-Bowie  through-crack 
solution  until  final  fracture. 

B.  Newman-Bowie  (New  model) 

This  program  is  similiar  to  the  program  in  A.  above. 
From  back  surface  penetration  until  a  through-crack 
develops  (i.e.,  c'/c=1),  the  Grandt-Bowie  solution  is 

corrected  with  the  transition  corrections  developed  in  this 
study.  These  factors  are  functions  of  c'/c.  After 


c'/c  =  0.9>  the  uncorrected  Grandt-Bowie  solution  for  a 
through-crack  is  used  until  final  fracture. 

C.  Newman-Bowie  (Opel's  model) 

This  program  uses  a  Newman- Ra  j  u  solution  until  a/t=0.75. 
At  that  point,  Opel  hole  correction  factors  are  applied  to 
the  Newman-Raju  solution.  The  Opel  model  uses  0  and  90  for 
the  parametic  angle,  .  After  back  surface  penetration, 
until  c/R=2.5»  surface  correction  factors  are  applied  to 
the  Grandt-Bowie  solution.  From  c/R=2.5  until  final 
fracture,  the  uncorrected  Grandt-Bowie  solution  is  used. 

D.  Newman-Bowie  ( Colli priest-Ehret  model) 

This  program  uses  an  uncorrected  Newman-Raju  corner- 
crack  solution  until  back  surface  penetration.  After  back 
surface  penetration  until  the  end  of  the  transition  region, 
an  imaginary  crack  length  (FSC)  is  used  in  the  Newman-Raju 
solution  to  calculate  the  front  surface  stress  intensity 
factors.  A  correction  factor  (CEC)  is  applied  to  the 
Grandt-Bowie  through-crack  solution  to  calculate  the  back 
surface  stress  intensity  factors.  The  transition  region 
ends  when  the  front  and  back  surface  crack  lengths  are 
equal.  After  this  point,  the  uncorrected  Grandt-Bowie 
solution  is  used  until  final  fracture. 

E.  Brussat 

This  program  uses  an  engineering  approximation  developed 
for  a  quarter-circular  crack-at- a-hol e.  It  uses  an  initial 
crack  length  that  is  the  average  of  the  actual  initial  test 


specimen  crack  lengths  along  the  bore  of  the  hole  and  along 
the  front  surface.  This  model  uses  the  Newman-Raju  finite 
width  correction. 

F.  Wal ker / Grandt- Bow  ie  Surface  Stress  Intensity 
Corr el  a  tions 

This  program  calculates  the  experimental  stress 
intensity  factors  based  on  curve-fit  data  for  c  vs  dc/dN 
experimental  data.  It  also  calculates  analytical  stress 
intensity  factors  using  the  Grandt-Bowie  through-crack 
solution.  The  front  and  back  surface  crack  length  ratio, 


PROGRAM  TRNSUN3(  INPUT,  OUTPUT,  TAPE  5  =  IN  PUT,  TAPE6  =  OU  TPU  T  ) 
INTEGER  X,XX 

REAL  PI ,  PA,  PC,  K9,  KA,  KC,  N,  NC,  LI  ,L2  ,M1  ,M2  ,M3  ,  A,  C,  R,  S,  T,  W, 

&  AS  TO  P , RS 
DO  40  1=1,6 
K9=53  .5 

PRINT*,  ’  INPUT  MATERIAL  AND  GEOMETRY  VARIABLES’ 

READ  (  5  ,  * )  A,  C,  R, S, T, W, ASTOP, RS 

IF(RS.GT.O.O)  GO  TO  5 

C1=4. 31544 E- 09 

N=3  -21 16 

FAC=0 . 0 

GO  TO  6 

5  C1  =  1  .  1  1616E-09 
N=3 .71797 

FAC  =  0 .7 

6  CONTINUE 

WRITE  ( 6 ,35)A, C, R, S, T, W, ASTOP, RS 
35  FORMAT ( 8 F 8 . 3 ) 

PA= 1 .39626 
PC=0 .17453 
NC=  1  .0 
B=W/2  . 

PI=3 . 1 41 6 
XX=1 
X  =  1 

WRITE ( 6 , 130) 

This  is  the  Newman-Raju  corner  crack  solution  with  the 
parametric  angle  =  1 0  &  80 

10  F 1  =  C/A 
F2= A/T 
Z3=A/C 

LI =( 1 ./(  1  .  +  ( C / R ) *COS (  .85* (PA) ) ) ) 

L2  =  (1  ./(I  . +( C/R ) *COS(  .85* (PC)))) 

Q 1  = (  1  .  +  1  . 46  4  *F 1 ** 1.65) 

M1  =  (F1**.5)*(1  .  +  .  0  4  *F  1  ) 

M2  =  ( . 2  *F 1 **  4  .  ) 

M3  =  (- . 1 1 *F1**4  .  ) 

G1=( 1 .+( .1+.35*F1*F2**2. )*(1 .-SIN(PA))**2. ) 

G2=( 1 . +( . 1+ . 35*F1 *F2**2 . )*( 1 .-SIN (PC ) )**2 . ) 
G3=((1.-.15*L1+3.46*L1*«2.-4.47*L1»*  3 • +3 . 52*L 1 *• 4 . ) / ( 1 . 

4  +  .  13*L1*»2 . ) ) 

G4=( ( 1 .- . 1 5*L2+3.46*L2**2.-4 .47*L2**3 .+3 -52*L2**4. )/( 1 . 

4+ . 1 3*L2**2 . ) ) 

G5=(1.13-.09#F1)*(1.+.1»(1.-C0S(PA))**2.)*(.8+.2*F2**.25) 
G6  =  ( 1  . 1 3- .0  9*F1 )•( 1  .  +  .  1*( 1  ,-C0S( PC ) ) **2 . )•(  . 8+. 2*F2**.25) 
F3  =  ( ( F 1  *  *  2 .  )*(SIN(PA))**2.+(C0S(PA) )**2. )**.25 
F4=((F1**2.)*(SIN(PC))*»2.+(C0S(PC))**2.)**.25 
F  5  =  (  1/ (COS( ( PI*R)/( 2. *B) ) ) ) 

F6  =  ( 1/ (COS( (PI*(2 . *R+NC«C )/( 4. »( B-C )+2 . *NC*C ) ) *F2*«.  5) ) ) 


F7=(F5*F6)«».5 

F8= (M 1 +M2*F2**2 ,+M3*F2**4. ) *G 1 *G3 *G5 *F3*F7 
F9=(M1+M2«F2«*2 .+M3*F2**4 . ) «G2»G4*G6*F4*F7 
D1=(S*SQRT(PI*A/Q1 ) )*F8 
D2=(S*SQRT( PI*A/Q1 ) ) »F9 

KA=D1*(((4./PI+(A*C)/(2.*T*R))/(4./PI+(A*C)/(T*R)))**.5) 
KC=D  2* ( ( (4 ./PI+(A*C)/(2.  »T*R ) ) / ( 4 . / PI+ ( A «C ) / ( T*R ) ) ) «• . 5 ) 
130  F0RMAT(2X, ‘A*  ,5X,  'C*  ,5X,  'A/T*  ,3X,  'A/C*  ,6X,  «KA«  ,6X, *KC’  ,9 
AX, » DA' , 9X, » DC* , 9X, ’CYCLES') 

DA=( ( (KA)*(1-RS)**(FAC))**N)*C1 
DC= ( ( (KC)*( 1-RS)«*(FAC) )**N)*C1 

IF ( (X/( XX*500) ) . EQ. 1 .0)  WRITE(6,15)  A , C , F2 , Z3 , KA , KC , 

ADA,  DC,  X 

IF((X/(XX*500)) .EQ.1 .0)  XX=XX+1 

IF(X.EQ.I)  WRITE (6, 15)  A , C , F2 , Z3 , KA , KC , DA , DC , X 

IF(KA.GE.K9)  GO  TO  200 

IFCKC.GE.K9)  GO  TO  200 

A= A+DA 

IF  (A.GT.T)  GO  TO  200 

This  signals  back  surface  penetration.  The  remainder 
of  the  program  is  an  uncorrected  Grandt-Bowie  through- 
crack  solution. 

C=  C+DC 
X  =  X  +  1 
GO  TO  10 

15  FORMAT  (3F6.3,F7.3,2F10.2,2E11.4,I8) 

00  WRITE (6, 15)  A, C,F2, Z3, KA, KC, DA, DC, X 
71  A=  C 

FWC=(COS(PI*(R+A)/4. ))••(-. 5) 

AA=(. 6762062 +(. 87 33015/C . 3 245 4 42+ A/ R ) ) ) »S» ( PI«A ) •• . 5 
KA=  AA  *FW  C 
KC  =  KA 

DA=( ( ( KA ) • ( 1-R3)**(FAC))**N)*C1 
DC=( ( ( KC ) * ( 1 - RS ) ** ( FAC ) )**N)*C1 

IF((X/(XX*500)).EQ.1 .0)  WRITE{6,15)  A , C , F2 , Z3 , KA , KC , 

ADA, DC, X 

IF  ((X/(XX*500)) .EQ.1 .0)  XX=XX+1 

IF  (X.EQ.1)  WRITE  (6,15)  A , C , F2 , Z3 , K A , KC , DA , DC , X 

IF(KA.GE.K9)  GO  TO  600 

IF  (KC.GE.K9)  GO  TO  600 

A=  A+  DA 

C=  C+DC 

X  =  X+1 

GO  TO  71 

600  WRITE  (6,15)  A  ,  C ,  F2 , Z3 , KA , KC , DA , DC  ,  X 
40  CONTINUE 


PROGRAM  TRAN  COR  (IN  PUT,  OUTPUT,  TAPE5  =  INPUT,  TA  PE6  =  OU  TPU  T  ) 
INTEGER  X,XX 

REAL  PI ,  PA,  PC,  K9,  KA,  KC,  N,  NO,  LI  ,  L2  ,M1  ,M2  ,M3  ,  A,  C,  R,  S,  T,  W, 

4  AS  TO  P , RS 
DO  40  1=1,6 
K9=53  .5 

PRINT*,  »  INPUT  MATERIAL  AND  GEOMETRY  VARIABLES’ 

RE  AD  (  5  ,  *  )  A,  C,  R,  S,  T,  W,  ASTOP,  RS 

IF(RS.GT.O.O)  GO  TO  5 

C  1=4 . 31544E-09 

N=3 .21  16 

FAC  =  0 .0 

GO  TO  6 

5  Cl =1 . 1 16 16E-09 
N=3. 71797 

FAC  =  0 .7 

6  CONTINUE 

WRITE  (6 ,35)A, C, R, S, T, W, ASTOP, RS 
35  FORMAT ( 8F8 . 3 ) 

PA= 1  .39626 
PC=0 . 17453 
NC=  1  .0 
B=  W/2  . 

PI=3 • 1416 
XX=1 
X  =  1 

WRITE ( 6 ,130) 

This  is  the  beginning  of  the  Newman-Raju  solution, 
uncorrected,  with  the  parametric  angle  =  1 0  4  80 

10  F  1  =C/A 
F2= A/T 
Z  3  =  A/C 

LI  =  (  1  ./(  1 .+(C/R)*COS( .85* (PA) )  ) ) 

L2  =  (1  ./(I  . + ( C/ R ) *COS(  .85* (PC) ))) 

Q1  =(  1  .  +  1  . 4  6  4  *F  1  **  1.65) 

M1  =  (F1**.5)*(1  .  +  . 0  4  *F 1 ) 

M2  =  (  . 2  *F 1 •*  4 . ) 

M3  =  ( - • 1 1 *F 1 **4  .  ) 

G1 =( 1 .+( .1 +.35*F1*F2**2. )*( 1 .-SIN (PA) )  •* 2 . ) 

G2=( 1 .+( . 1+.35*F1*F2**2. )  *  (  1 .-SIN(PC) )**2. ) 

G3=( ( 1 .- . 1 5*L1 +3 .46*L1 *«2.-4 .47*L1 »*3.+3 .52#L1 »»4. ) / ( 1 . 
4+. 1 3*L1**2 . ) ) 

G4  =  ( ( 1 . - . 1 5*L2  +3 .46*L2**2.-4.47*L2**3.  +  3.52*L2**4. ) /(  1 . 

4+ . 1 3*L2«*2 . ) ) 

G5=(1.13-.09#F1)#(1.+.1*(1.-C0S(PA))**2.)»(.8+.2»F2»».25) 
G6  =  ( 1  . 1 3- .0  9#F  1  )•( 1  .  +  .  1*( 1  .-C0S( PC ) ) **2 .)*(.8+.2«F2«*.25) 
F3  =  ( ( F 1 **  2 . ) * ( SIN ( PA ) ) •• 2 .  +  ( COS ( PA ))»«2. )*»-25 
F4  =  ( ( F  1  *  *  2 .  )*(SIN(PC) )**2.+(C0S(PC) )**2. )**.25 
F5=( 1/(COS( (PI*R)/( 2. *B) ) ) ) 


F6=(1/(C0S((PI*(2.*R+NC*C)/(4.*(B-C)+2.*NC*C))*F2**.5))) 

F7=(F5*F6)*«.5 

F8=(M1+M2«F2**2.+M3#F2»»4.)»G1*G3*G5#F3*F7 
F9=(M1+M2*F2i#2.+M3*F2**4 . )*G2*G4*G6*F4*F7 
D1=(S*SQRT(PI*A/Q1 ) ) *F8 
D2=(S*SQRT(PI*A/Q1 ) )  *F 9 

K  A  =  D  1  *  ( ( ( 4./PI+(A*C )/( 2. *T*R ) )/( 4./PI+(A*C)/(T*R) ) )**. 5) 
KCsD2*(((4./PI+(A*C)/(2.*T*R))/(4./PI+(A*C)/(T*R)))**.5) 
130  FORMA T ( 2X ,  'A',5X,'C',5X,'A/T',3X,'A/C',6X,'KA',6X,'KC',9 
&X, ' DA' , 9X, • DC' , 9X, • CYCLES' ) 

DA=( ( ( KA ) * ( 1-RS)**(FAC) )**N )*C 1 
DC=( ( ( KC ) * ( 1-RS)**(FAC) ) *  *  N ) *C 1 

IF( (X/(XX»500) ) . EQ. 1 .0)  WRITE(6,15)  A  ,  C  ,  F2 , Z3 , KA , KC  , 

&DA , DC, X 

IF((X/(XX*500)).EQ.1 .0)  XX=XX+1 

IF(X.EQ.I)  WRITE (6,15)  A , C , F2 , Z3 , KA , KC , DA , DC , X 

IF(KA.GE.K9)  GO  TO  200 

IF ( KC . G E. K 9 )  GO  TO  200 

A= A+DA 

IF  (A.GT.T)  GO  TO  200 

This  signals  back  surface  penetration  and  the  Chansler 
correction  factors  are  applied  to  the  Grandt-Bowie 
through-crack  solution  until  the  end  of  the  transition 
region 

C=  C+DC 
X  =  X+1 
GO  TO  10 

15  FORMAT  (3F6.3,F7.3,2F10.2,2E11 .4,18) 

30  WRITE (6, 15)  A, C, F2 , Z3 , KA, KC, DA, DC, X 
A=  .005 
C=  C+  A 

61  FWCs(C0S(PI»(R+A)/4. ))••(-. 5) 

AAs( .6762062+(.8733015/( . 3 245 442+ A/ R ) ) ) »S» ( PI *A ) •• . 5 

BCF1 =- .6896 1 4* A/C 

BCF2  =  1  . 2327  5* ( A/C ) **2 

BCF= . 609633+BCF1+BCF2 

KA=  AA  *FWC/BCF 

FW C=(C OS (PI*(R+C)/4. ))**(-. 5) 

AA=(.6762062+(.8733015/(.3245442+C/R)))«S»(PI#C)**.5 
SCF 1 =- . 1 3564*A/C 
SCF  =  1  .  24007  +  SCF 1 
KC  =  AA*FWC/ SCF 

DA=(  (  (  K  A  )  •  (  1-RS)**(FAC))**.N)*C1 
DCs( ( (KC)»(1-RS )*•( FAC ))••») «C1 

IF( (X/(XX*500) ) . EQ. 1 .0)  WRITE(6,15)  A , C , F2 , Z3 , KA , KC , 

&DA, DC, X 

IF  ( (X/(XX*500) ) . EQ. 1 .0)  XX=XX+1 

IF  (X.EQ.1)  WRITE  (  6  ,  1  5  )■  A  ,  C  ,  F2  ,  Z3  ,  KA,  KC,  DA,  DC,  X 
IF  (KC. GE. K9)  GO  TO  600 
A=  A+DA 


C=  C+  DC 

IF.  (  (A/C)  .GE.O  .9)  GO  TO  400 

This  signals  the  end  of  the  transition  region 
The  remainder  of  the  program  is  an  uncorrected 
Grandt-Bowie  through-crack  solution 


X  =  X  +  1 
GO  TO  61 

400  WRITE  (6,15)  A , C , F2 , Z3 , KA , KC , DA , DC , X 
71  A=C 

FWC= ( COS( PI*(R  +  A)/4 . ))••(-.  5) 

A A=  (  .676 206 2  +  ( .873  301  5/ ( . 3 24 5 4 42  + A/ R ) ) ) *S« ( PI *A ) ** . 5 
KA= AA*FW C 
KC  =  KA 

DA=( ( ( KA ) * ( 1-RS)**(FAC))**N)*C1 
DC=( ( ( KC ) * ( 1-RS)**(FAC) ) **N ) *C 1 

IF( (X/(XX«500) ) . EQ. 1 .0)  WRITE(6,15)  A , C , F2 , Z3 , KA , KC , 
&DA, DC, X 

IF  ( (X/(XX«500) ) .EQ.1 .0)  XX=XX+1 

IF  (X.EQ.1)  WRITE  (6,15)  A , C , F2 , Z3 , KA , KC , DA, DC , X 

IF(KA.GE.K9)  GO  TO  600 

IF  (KC.GE.K9)  GO  TO  600 

A=  A+DA 

C=  C+  DC 

X  =  X+1 

GO  TO  71 

600  WRITE  (6,15)  A , C , F2 , Z3 , KA , KC , DA , DC , X 
40  CONTINUE 
END 


PROGRAM  TR  AN  S3  (  INPUT,  OUTPUT,  TAPE  5  =  IN  PUT,  TAPE6  =  OU  TPU  T ) 
INTEGER  X,XX 

REAL  PI,  PA,  PC,  K9.KA,  KC,  N,  NC,  LI  ,L2  ,M1  ,M2,M3  ,  A,  C,  R,  S,  T,  W 
&  AS  TO  P , RS 
DO  40  1=1,6 
K9=53  .5 

PRINT*,  '  INPUT  MATERIAL  AND  GEOMETRY  VARIABLES' 

READ ( 5 , *  )  A,  C,  R,  S, T, W, ASTOP, RS 

IF(RS.GT.O.O)  GO  TO  5 

Cl =4 . 31  544E-09 

N  =  3 .21  16 

FAC  =  0  .0 

GO  TO  6 

5  Cl  =  1  .  1  16  16E-09 
N=3. 71797 

FAC  =  0 .7 

6  CONTINUE 

WRITE  (6 ,35) A, C, R, S, T, W, ASTOP, RS 
35  FORMAT ( 8F8 . 3 ) 

PA=  1  .5707  96 
PC  =  0  .0 
NC=  1  .0 
B=W/2  . 

PI=3 • 1 41 6 
XX=  1 
X  =  1 

WRITE(6,130) 

This  begins  the  Newman-Raju  corner-crack  solution 

10  FI =  C / A 
F2= A/T 
Z3=A/C 

LI =(  1  ./( 1 .+(C/R)*C0S( .85* (PA) ) ) ) 

L2=(1 ./(I .+(C/R)*COS( .85* (PC) ) )) 

Q1  =  (  1  .  +  1  .464*F1**1 .65) 

Ml  = ( F 1 ** . 5  )  * (  1  .  +  . 0 4 *F  1  ) 

M2  =  (  . 2  *F 1 **  4  .  ) 

M3  =  ( - • 1 1 *F 1 **4  .  ) 

G1=(1.+(.1+.35*F1*F2*«2.)*(1 .-SIN( PA ) ) *»2. ) 

G2=( 1 . +( . 1+. 35*F1 *F2**2 . )*( 1 .-SIN( PC) ) **2 . ) 
G3=((1.-.15*L1+3.46*L1**2.-4.47*L1**3.+3.52*L1**4.)/(1 
4+. 1 3  *L 1 **2 . 1)) 

G4  =  ( ( 1  .  - . 1 5*L2  +  3 .46*L2**2 .-4 .47*L2**3 .  +  3 .5  2*L2**4 . ) /(  1 
4  + . 1 3*L2*«2 . 1 ) ) 

G  5  =  ( 1 • 1 3 - • 0  9*F 1 ) *(1.  +  .1*(1.-C0S(PA))**2.)*(.8  +  .2*F2«* 
4.25) 

G6=(1.13-.09*F1)*(1.+.1*(1.-COS(PC))**2.)*(.8+.2*F2*« 

4.25) 

F3=((F1**2.)*(SIN(PA))»«2.+(C0S(PA))**2.)**.25 

F4=((F1**2.)*(SIN(PC))*»2.+(C0S(PC))**2.)**.25 


F  5  =  (  1/(C0S( ( PI »R ) / ( 2 . *B ) ) ) ) 

F6=( 1/ (COS( ( PI*( 2. »R+NC«C )/( 4. *(B-C)+2. «NC»C ) ) »F2»«. 5) ) ) 
F7=(F5*F6)**.5 

F8=(M1+M2«F2»*2. +M3#F2»»4. ) »G 1 *G3 «G5 *F3*F7 
F9=(M1+M2*F2*»2 .+M3#F2*»4 . ) *G2»G4«G6»F4*F7 
D1=(S«SQRT(PI*A/Q1 ) )»F8 
D2=(S«SQRT(PI»A/Q1 ) ) *F9 
IF  ( F2 . LE. . 75 )  GO  TO  100 

If  this  conditon  is  met,  Opel's  hole  bore  corrections  take 
effect  here 

TF 1  =  1  . 4  6  80  6  8*F  2 
TF2  =  - 1  . 27  57  87  *F2**2 
TCF= .573924+TF1+TF2 

KA  =  D 1  * ( ( ( 4./PI+(A*C)/( 2. *T*R) ) /  (  4 . / PI+ ( A «C ) / ( T»R ) ) ) •*  .  5  ) 
&/TCF 

SCI =- .362863*F2 
SCF  =  1  .3  16  867  +  SC 1 

KC  =  D2*( ( ( 4./PI+( A*C)/(2. #T«R) )/( 4 . / PI+ ( A »C ) / ( T«R )  ) )•». 5) 
&/SCF 
GO  TO  101 

100  KA  =  D1*(((4./PI+(A»C)/(2.*T*R))/(4./PI+(A»C)/(T«R  )))*«. 5) 

KC  =  D  2* ( ( (4 ,/PI+(A*C) /(2. *T*R) )/ (4 ./PI+(A*C)/ ( T  *R ) ) ) **. 5) 

101  CONTINUE 

130  FORMAT (2X,  *  A »  , 5X ,  » C»  , 5X,  • A/T *  , 3X,  • A/C'  , 6 X ,  • KA '  , 9X ,  ' KC  '  , 9 
&X,  ' DA'  , 9X,  '  DC'  , 9X, ' CYCLES' ) 

DA=( ( ( K A  )  • ( 1-RS)**(FAC))**N)*C1 
DC=( ( (KC)*( 1-RS)**(FAC) )**N)*C1 

IF((X/(XX»500) ) .EQ.1 .0)  WRITE(6,15)  A , C , F2 , Z3 , KA , KC , 

&DA,  DC,  X 

IF((X/(XX*500)) .EQ.1 .0)  XX=XX+1 

IF(X.EQ.I)  WR ITE ( 6 , 1 5 )  A , C , F2 , Z3 , K A , KC , DA , DC , X 

I F ( KA . GE. K9)  GO  TO  200 

IF ( KC. GE. K9)  GO  TO  200 

A=  A+  DA 

IF  (A . GT. T)  GO  TO  200 

If  this  condition  is  met,  back  surface  penetration  has 
occurred  and  Opel's  surface  corrections  are  applied  to 
the  Grandt-Bowie  through-crack  solution 

C=  C+  DC 
X  =  X  +  1 
GO  TO  10 

15  FORMAT  ( 3F6.3 ,F7.3 ,2F10.2,2E1 1 .4 , 18) 

00  WRITE (6, 15)  A, C, F2 , Z3 , KA, KC, DA, DC, X 
61  A=R*(-1 .931 1+1 .70076*C/R) 

FWC=(COS(PI*(R+A)/4 .))«*(  -.5) 

AA=(.676  20  62  +  (.  87  33015/(.3245442-t-A/R)))*S»(PI»A)«».5 

BCF1 = .72496 1 *A/R 

BCF2  =  - . 1  58451  * ( A/R ) **2 


BCF= . 1 916C1+BCF1+BCF2 
KA=  AA«FWC/BCF 

FWC=(C0S(PI»(R+C)/4. ))••(-. 5) 

AA=(.6762062+(.8733015/(.3245442+C/R)))*S»(PI»C)*«.5 

SCF1 =- . 1 06745 *C/R 

SCF2= . 0  1  8451  * (C/R ) *»2 

SCF= 1  .27  9537+SCF1+SCF2 

KC  = AA*FWC/SCF 

DA=(  (  ( KA ) * ( 1-RS)**(FAC) ) **N ) *C  1 
DC= ( ( ( KC ) * ( 1 - RS ) «• ( FAC ) ) **N ) «C 1 

IF( (X/(XX*500) ) . EQ. 1 .0)  WRITE(6,15)  A , C , F2 , Z3 , KA , KC , 
&DA,  DC,  X 

IF  ( (X/(XX«500) ) . EQ. 1 .0)  XX=XX+1 

IF  (X.EQ.1)  WRITE  (6,15)  A , C , F2 , Z3 , K A , KC , DA , DC , X 
IF  (KC.GE.K9)  GO  TO  600 
A=  A+DA 
C=  C+  DC 

IF  ( (A/R) .GT.2.5)  GO  TO  400 

Meeting  this  condition  signals  the  end  of  Opel's 
transition  region  and  the  rest  of  the  program 
is  an  uncorrected  Grandt-Bowie  through-crack 
solution 

X  =  X+1 
GO  TO  61 

400  WRITE  (6,15)  A  ,  C ,  F2  ,  Z3  ,  KA  ,  KC  ,  DA  ,  DC  ,  X 
71  A=  C 

FWC=(C0S(PI*(R+A)/4. ))•*(-. 5) 

AA=(  .6762062  +  (  .  873301  5/ (  .  32 45 4 42  + A/R ) ) ) »S* ( PI »A ) •• . 5 
KA=  AA*FW  C 
KC  =  KA 

D  A  =  ( ( ( K A  )  * ( 1-RS)**(FAC) ) *  *N ) *C  1 
DC=( ( ( KC ) * ( 1  - RS ) ** ( FAC ) )**N)«C1 

IF  ( (X/(XX*500) ) .EQ. 1 .0)  WRITE  (6,15)  A , C , F2 , Z3 , KA , KC 
&DA , DC, X 

IF  ( (X/(XX*500)) .EQ. 1 .0)  XX=XX+1 

IF  (X.EQ.1)  WRITE  (6,15)  A , C , F2 , Z3 , K A , KC , DA , DC , X 

IF  (KA.GE.K9)  GO  TO  600 

IF  ( KC. GE. K9)  GO  TO  600 

A=  A+DA 

C=  C+DC 

X  =  X  +  1 

GO  TO  71 

600  WRITE  (6,15)  A , C , F2 , Z3 , KA , KC , DA , DC  ,  X 
40  CONTINUE 


PROGRAM  COL  EH  R2 ( I N PU T , OU TPU T ,  TA PE 5  =  IN PU T ,  TAPE 6  =  OU TPU T ) 
INTEGER  X,XX 

REAL  PI ,  PA  ,  PC,  K9  ,  KA  ,  KC  ,  N,  NC,  LI  ,  L2  ,  Ml  ,  M2  ,  M3  ,  A,  C  ,  R,  S,  T,  W 
&  AS  TO  P ,  RS 
DO  40  1=1,6 
K 9=  53  .5 

PRINT*,  '  INPUT  MATERIAL  AND  GEOMETRY  VARIABLES’ 

READ  (  5  ,  *  )  A,  C,  R,  S,  T,  W,  ASTOP,  RS 

IF( RS. GT. 0 .0 )  GO  TO  5 

Cl =4 . 3  1  544  E-0  9 

N=3 .2116 

FAC  =  0  .  0 

GO  TO  6 

5  Cl =1  . 1 1 6 1 6  E-0  9 
N=3 .7 1797 

FAC  =  0 .7 

6  CONTINUE 

WRITE  (6 ,35) A, C, R, S, T, W, ASTOP, RS 
35  FORMAT(  8F8 .3) 

PA=  1  .39626 
PC=0 .17453 
NC= 1 .0 
B=W/2 . 

PI  =  3 • 1 41  6 
XX=  1 
X  =  1 

WRITE ( 6 ,130) 

This  is  the  beginning  of  the  uncorrected  Newman-Raju 
corner-crack  solution  with  parametric  angle  =  10  &  80 

10  F 1 =  C/ A 
F2= A/T 
Z3  =  A/C 

L1  =  (  1  .  / (  1.  +  (C/R)*COS(  .85*(PA)))) 

L2  =  ( 1  .  /( 1 . +(C/R )*C0S( . 85* (PC) ) ) ) 

Q1  =  (  1  .  +  1  . 4  6  4  *F 1  *  *  1 .65) 

Ml  =  (F 1 *• . 5 ) *( 1  .  +  .  0  4  *F  1  ) 

M2  =  (  . 2  *F  1  **  4  .  ) 

M3  =  ( - . 1  1  *F  1  **4  .  ) 

G1 =(  1  .  +  (  . 1  +  .35*F1*F2**2. )*(  1  .-SIN( PA) )  *  *  2  . ) 

G2  =  ( 1  .+(  . 1  +  . 3  5  #F 1 *F 2 * *2  . )  *  (  1  .-SIN(PC) ) **2 . ) 
G3=((1.-.15*L1+3.46*L1«»2.-4.47«L1**3.+3-52*L1*»4.)/(1 
&+ . 1 3  *  L 1  * *2 . ) ) 

G4= ( ( 1 .- . 1 5*L2+3 .46*L2*«2.-4 .47*L2»*3 .+3 .52*L2**4 . ) /( 1 
&+ . 1 3*L2**2 . ) ) 

G  5  =  (  1  •  1 3  —  -  09*F1 ) * ( 1  .  +  . 1  * ( 1 .-COS( PA ) )**2. ) •(  .8+.2*F2** 

&  •  2  5  ) 

G6=(1.13-.09*F1)«(1.+.1»(1.-COS(PC))»*2.)*(.8+.2»F2*» 

4.25) 

F  3  =  (  C  F 1  *  *  2 .  ) * ( S  IN ( PA) )**2.  +  (C0S( PA) ) • *  2 . )**.25 


F  4  =  (  ( F  1  *  *  2 . ) * ( S IN ( PC) )  **2 .  +  ( COS (PC ) ) »*2 . ) ••  .25 
F  5  =  (  1  / ( C  0  S ( (  P I *R  )  /  ( 2 . *B  )  ) ) ) 

F6=( 1 / ( COS ( (PI*(2.*R+NC*C)/(4.*(B-C)+2.*NC*C))*F2#*.5)) ) 
F7=(F5*F6)««.5 

F8=(M1+M2»F2«*2.+M3#F2**4. )«G1 «G3#G5*F3*F7 
F9= (M1+M2*F2**2 .+M3*F2*«4 . ) *G2 »G4 »G6 »F 4 *F7 
D1=(S*SQRT(PI*A/Q1 )  )  *F8 
D2=(S«SQRT( PI»A/Q1 ) ) «F9 

KA=Dt*(((4./PI+(A*C)/(2.*T*R))/(4./PI+(A*C)/(T*R)))**.5) 
KC=D2*(((4./PI+(A*C)/(2.*T#R))/(4./PI+(A*C)/(T*R)))**.5) 
130  FORMA T(2X,  *  A '  , 5X, • C ' , 5X ,  ' A/T • , 3 X,  »A/C'  , 6 X ,  • KA '  , 6 X ,  ' KC '  , 9 
&X, ' DA' , 9X, ' DC' ,9X, 'CYCLES') 

DA=( ( ( KA ) * ( 1 - RS ) »• (FAC ) ) **N ) *C 1 
DC= ( ( (KC)*(1-RS)**(FAC))**N)*C1 

IF( (X/(XX*500) ) . EQ. 1 .0)  WRITE(6,15)  A , C , F2 , Z3 , KA , KC , 

&DA, DC, X 

IF( (X/(XX*500) ) .EQ. 1 .0)  XX=XX+1 

IF(X.EQ.I)  WRITE (6, 15)  A, C , F2 , Z3 , K A , KC , DA , DC , X 

IF( KA . GE. K9)  GO  TO  200 

IF  (  KC.  GE.  K9)  GO  TO  200 

A= A+DA 

IF  (A.GT.T)  GO  TO  200 

This  signals  back  surface  penetration.  The  Collipriest- 
Ehret  corrections  are  applied  from  here  to  the  end  of  the 
transition  region.  An  imaginary  crack  length  (FSC)  is 
used  with  the  Newman-Raju  corner-crack  solution  to  calcu¬ 
late  the  front  surface  stress  intensity  factors  and  a 
correction  factor  ( CEC )  is  applied  to  the  Grandt-Bowie 
through-crack  solution  to  calculate  the  back  surface 
stress  intensity  factors. 


C= C+DC 
X  =  X+1 
GO  TO  10 

15  FORMAT  (3F6.3,F7.3,2F10.2,2E11 .4,18) 

DO  WRITE (6, 15)  A, C , F2 , Z3 , KA, KC , DA , DC , X 
CB= .020 
CF  =  C 

P A  =  0 .17453 

7 1  FSC=  T* ( ( 1 .0-(CB/CF)«* 2.0 )«•(-. 5) ) 

A=  FSC 

FWC=(C OS (PI*(R+CB)/4. ))**(-. 5) 

AA=(  .676  2062  +  (  .  87 3  3  0 1  5 / (  .  3245442  +  CB/R) ) ) *S • ( PI »CB ) * * . 5 
CEC=( 1 . 0 / (  1  .0-( 1  .0-(CB**2.0/CFt#2.0) ) **.5) ) •• . 5 
FI =CF/ A 
F2= 1 .0 
Z3= A/CF 

L1=(1./(1.+(CF/R)*COS(.85#(PA)))) 

Q1  =  (  1  . +1  .  46  4  #F 1 ** 1 .65) 

M1  =  (F1**.5)*(1 •  +  . 0  4*F 1 ) 


M2  =  (  .  2  *F 1 **4  .  ) 

M3  =  ( - • 1 1 *F 1  *  *  4  .  ) 

G1 =(  1  .  +  (  .1  +  .35*F1*F2**2. )•( 1 .-SIN (PA) )**2 . ) 

G3  =  ( ( 1  • 1 5*L1  +3 . 46»L1 *»2 .-4 . 47 *L 1**3 .+3 . 52 *L 1**4 . )/( 1 

&+. 13*L1«*2. ) ) 

G5  =  ( 1  .1 3-.0  9*F1 ) •( 1 .  +  . 1 *( 1 .-COS( PA) ) •* 2 . ) * ( . 8+ . 2 *F2** 

&  .  2  5  ) 

F  3  =  C (F1**2.)*( SIN (PA))**2.+(C0S(PA))**2.)**.25 
F5  =  (  1/(C0S((PI»R)/(2.*B)))) 

F6  =  ( 1 / ( COS ( ( PI  * ( 2 . *R+NC*CF)/( 4. *(B-CF)+2.*NC*CF) ) *F2*» 
4.5)  )  ) 

F7=(F5*F6) *». 5 

F8=(M1+M2*F2**2.+M3«F2**4. )»G1 *G3*G5*F3*F7 
D1 =  ( S  *SQRT ( PI  *A  /  Q 1 ) ) #F8 

KC=D1*( ( (4./PI+(A#CF)/(2.*T*R) )/( 4./PI+(A*CF)/(T*R) ) ) 

4  **  .  5 ) 

KA= AA*FWC*CEC 

DA=  (  (  (KA)*(1-RS)**(FAC))**N)*C1 
DC= ( ( ( KC ) * ( 1 - RS ) • • ( FAC ) )»»N)»C1 

IF((X/(XX*500)) .EQ. 1 .0)  WRITE(6,15)  A , CF , F2 , Z3 , KA , KC , 
&DA,  DC,  X 

IF  ( (X/(XX»500) ) . EQ. 1 .0)  XX=XX+1 

IF  (X.EQ.1)  WRITE  (6,15)  A, CF, F2 , Z3 , KA, KC, DA, DC, X 
CB=  CB  +  DA 
CF  =  CF+  DC 
X  =  X  +  1 

IF ( CB . GE. CF)  GO  TO  400 

When  this  conditon  is  met,  the  front  and  back  surface 
crack  lengths  are  equal  and  this  signals  the  end  of 
the  transition  region.  The  rest  of  the  program  is  an 
uncorrected  Grandt-Bowie  through-crack  solution. 

GO  TO  71 

400  WRITE  (6,15)  A  ,  CF ,  F2  ,  Z3  ,  KA  ,  KC  ,  DA  ,  DC  ,  X 
75  A=  CF 

FWC=(COS(PI*(R+A)/4. ) )**(-. 5 ) 

AA=( .6762062+( . 873301 5/ ( . 3245442+A/R ) ) ) ( PI«A) . 5 
KA=  AA*FW  C 
KC  s  KA 

DA=( ( ( K A  )  • ( 1-RS)**(FAC))**N)*C1 
DC= ( ( (KC)»( 1-RS)«»(FAC) )*«N)»C1 

IF  ( (X/(XX*500) ) . EQ. 1 .0)  WRITE  (6,15)  A , CF , F2 , Z3 , KA , KC 
&DA, DC, X 

IF  ( (X/(XX*500) ) . EQ. 1 .0)  XX=XX+1 

IF  (X.EQ.1)  WRITE  (6,15)  A, CF, F2 , Z3 , KA, KC, DA, DC, X 

IF  (KA. GE. K9)  GO  TO  600 

IF  (KC.GE.K9)  GO  TO  600 

A=  A+DA 

CF=  CF+  DC 

X  =  X  +  1 


GO  TO  75 

600  WRITE  (6,15)  A , OF , F2 , Z3 , KA , KC , DA , DC , X 
40  CONTINUE 
END 


PROGRAM  B RUSS ( INPUT,  OUTPUT,  TAPE5  =  INPUT,  TAPE6  =  OU TPU T ) 
INTEGER  X,XX 

REAL  PI,  PA,  PC,  K9,  KA,  KC,  N,  NC,  LI  ,L2  ,M1  ,M2  ,M3,  A,  C,  R,  S,  T,  W 
&  AS  TO  P , RS 
DO  40  1=1  ,6 
K9=53 .5 

READ  (5,*)  A, C, R, S, T, W, ASTOP, RS 

IF  (RS.GT.O .0)  GO  TO  5 

C1=4.31544E-09 

N=3 .21  16 

FAC=0 .0 

GO  TO  6 

5  C 1  =  1  -  11616 E- 09 
N  =  3 .71797 
FAC=0 .7 

6  CONTINUE 

WRITE  (6,35)  A, C, R, S, T, W, ASTOP, RS 
35  FORMAT  (8F8.3) 

PI=3 • 1 416 

XX=1 

X=1 

WRITE  (6,130) 

A= ( C+ A  )  / 2 
C=  A 

10  F 1 =  C / A 
F2= A/T 
Z3=A/C 

FAC1=S*(PI*A)**.5 

FAC2= ( 1/C0S(PI*(A+2»R)/(2*(W-A) )))»•. 5 
FAC3  =  1-(  . 2  886/ (  1+2*(C/T)**2) ) 

FAC4=1. 2133-2. 205*(A/(A+R))+.6451*(A/(A+R))**2 
FAC5  =  E  XP ( FAC4 ) 

KA=FAC1*FAC2*FAC3*FAC5 
KC  =  KA 

30  FORMAT  ( 2X ,  ' A '  , 5 X ,  ' C '  , 5 X , ' A/T '  , 3X ,  ' A/C  *  , 6 X ,  ' KA '  , 6 X ,  ' KC 
&9X, ’ DA' , 9X, • DC* , 9X, ’CYCLES' ) 

DA=( ( ( KA ) * ( 1-RS)**(FAC) ) *  *N ) *C  1 
DC=( ( (KC)»(1-RS)»»(FAC))«*N)»C1 

IF  ( (X/( XX«500) ) .EQ.1 .0)  WRITE(6,15)  A , C , F2 , Z3 , KA , KC , 
4DA, DC, X 

IF  ( ( X/ ( XX*500 ) ) . EQ. 1 .0)  XX=XX+1 

IF  (X.EQ.1)  WRITE ( 6 , 1 5 )  A , C , F2 , Z3 , K A ,  KC ,  DA ,  DC ,  X 

IF( KA . GE.  K9)  GO  TO  200 

IF(KC.GE.K9)  GO  TO  200 

A=  A+DA 

C=  C+  DC 

X  =  X  +  1 

GO  TO  10 

15  FORMAT(3F6.3,F7.3,2F10.2,2E11 .4,18) 

200  WRITE ( 6 , 1 5 )  A , C , F2 , Z3 , K A , KC , D A , DC , X 
40  CONTINUE 

e: 


PROGRAM  INT7T(  INPUT,  OUTPUT,  TAPE5=  INPUT,  TAPE6  =  OU  TPU  T ) 
REAL  PI,  N,  M,  KE,  KA,  R,  T1  ,  T2,C0  ,C1  ,C2,C3,  C4,  C5  ,  B1  ,  B2,  MAXS 
PI=3 • 1 41 6 
DO  20  1=1  ,7 

PRINT*,  'INPUT  R,  COEFS  AND  MAX  STRESS' 

RE  AD  (  5  ,  *  )  R,  T1  ,  T2  ,  CO  ,  Cl  ,  C2  ,  C3  ,  C4  ,  C5  ,  B1  ,  B2  ,  MAXS 

IF  (R.GT.0.0)  GO  TO  100 

C=1  . 2 9  E- 0  8 

N  =  2 .89 

FAC  =  0 .0 

GO  TO  101 

100  C=3 . 2624 E- 09 
N=3 .3908 

FAC  =  0 .5 

101  CONTINUE 
D=  .20 

WRITE  (6,32) 

32  FORMAT  (4X, 'C' ,8X, 'C/C' ,8X, 'K  E XPER IMEN TAL * , 8X , 

4 ' K  ANAL YT' , 8X, » ACCURACY' ) 

DO  10  J=1 ,70 
A= . 0 1 5*  J+D 
T=T1  «A+T2 

This  part  uses  fifth-order  polynomial  regression  curve 
fit  data  based  on  experimental  results  to  calculate 
the  experimental  stress  intensity  factor  using  a  Walker 
equation. 

DADN=C0+C1*T+C2*T**2+C3*T»*3+C4*T**4+C5*T**5 

FACTOR=DADN/C 

COEF=ABS( FACTOR )•*( 1/N) 

KE=  COEF* ( 1 -R ) ** ( -FAC ) 

33  FORMAT ( 2X, F5. 3 ,5X, F7 . 3 , 5X , F8 . 3 , 1 5X , F8 . 3 , 1 2X , F6 . 3 ) 

This  part  calculates  the  analytical  stress  intensity 
factor  using  the  Grandt-Bowie  through-crack  solution 


RAD= .125 

FWC=(COS(PI*(RAD+A)/4. ))*•(-. 5) 

AA=  (  . 676 2062  +  (  .  87  3301  5/ ( . 32 454 42  + A/ RAD ) ) ) *MAXS* ( PI*A) 

4** .  5 
KA=  AA*FW  C 
ACCUR=KA/KE 

The  next  line  calculates  the  back  surface  crack  length 
based  on  linear  curve  fit  data  using  experimental  results 

B=B1  +B2  *  A 
AT=B/A 

10  WRITE(6 ,33)A, AT, KE, KA, ACCUR 
20  CONTINUE 
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This  study  develops  correction  factors  for  currently  used  stress 
intensity  factor  equations  to  more  accurately  predict  stress  intensity 
factors  for  a  corner-crack  emanating  from  a  hole  as  it  transitions  to 
a  uniform  through-the-thickness  crack.  These  correction  factors 
resulted  in  an  increase  in  accuracy  for  total  life  prediction  and  much 
better  correlation  between  analytical  stress  intensity  factor 
predictions  and  experimental  results  in  the  transition  region  for 
2024-T351  and  7075-T651  aluminum  alloys.  An  experimental  program  was 
undertaken  to  generate  all  2024-T351  aluminum  test  data, used  in  this 
investigation.  The  7075-T651  aluminum  test  data  was  generated  during 
earlier  work.  Correction  factors  developed  by  Opel  for  7075-T651 
aluminum  were  evaluated  for  2024-T351  aluminum  and  found  to  be 
excessively  conservative. 

Hartranft  and  Sih  suggested  the  stress  intensity  factors  be 
evaluated  at  an  imaginary  surface  at  an  angle  ty  away  from  the  front 
and  hole-bore  surfaces.  These  stress  intensity  factors  could  then  be 
used  in  the  life  prediction  models  for  life  from  a  corner-crack  until 
back  surface  penetration.  Evaluating  stress  intensity  factors  at  <J>= 
10°and  80°  eliminate  surface  boundary  effects  caused  by  fabrication 
processes  like  cold  rolling  and  hole  drilling,  thereby  improving  life 
predictions  to  back  surface  penetration  significantly.  It  is  shown 
that  with({>=  10°  and  80°  ,  no  other  corrections  are  necessary  for  life 
prediction  to  back  surface  penetration. 

A  transition  region  is  postulated  to  effect  life  prediction  from 
back  surface  penetration  until  final  fracture.  Correlations  were  made 
between  experimental  results  for  2024-T351  and  7075-T651  aluminum  test 
data  and  the  Grandt  linearization  of  the  Bowie  solution  for  a 
thorough-the-thickness  crack.  These  correlations  were  plotted  from 
back  surface  penetration  until  final  fracture.  The  plot  yields  the 
end  of  the  transition  region  and  the  required  correction  factors.  The 
life  prediction  model  incorporating  these  transition  correction 
factors  are  shown  to  be  the  most  accurate  and  versatile  of  all  models 
investigated. 

Life  predictions  were  made  using  the  new  model,  an  instantaneous 
model  (which  assumes  a  through-crack  at  back  surface  penetration), 
Opel's  model,  Brussat's  model  and  the  Collipriest-Ehret  model. 
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